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ABSTRACT 

We present subarcsecond resolution mid-infrared (mid-IR) photometry in the wavelength range 
from 8 to 20 ^m of eighteen Seyfert galaxies, reporting high spatial resolution nuclear fluxes for the 
entire sample. We construct spectral energy distributions (SEDs) that the AGN dominates, relatively 
uncontaminated by starlight, adding near-IR measurements from the literature at similar angular 
resolution. We find that the IR SEDs of intermediate-type Seyferts are flatter and present higher 10 
to 18 /^m ratios than those of Seyfert 2 galaxies. We fit the individual SEDs with clumpy dusty torus 
models using the in-house-developed BayesClumpy tool. We find that the clumpy models reproduce 
the high spatial resolution measurements. Regardless of the Seyfert type, even with high spatial 
resolution data, near- to mid-IR SED fitting poorly constrains the radial extent of the torus. For 
the Seyfert 2 galaxies, we find that edge-on geometries are more probable than face-on views, with a 
number of clouds along equatorial rays of A^o = 5-15. The 10 /im silicate feature is generally modeled 
in shallow absorption. For the intermediate-type Seyferts, A'o and the inclination angle of the torus are 
lower than those of the Seyfert 2 nuclei, with the silicate feature appearing in weak emission or absent. 
The columns of material responsible for the X-ray absorption are larger than those inferred from the 
model fits for most of the galaxies, which is consistent with X-ray absorbing gas being located within 
the dust sublimation radius whereas the mid-IR fiux arises from an area farther from the accretion 
disc. The fits yield both the bolometric luminosity of the intrinsic AGN and the torus integrated 
luminosity, from which we derive the reprocessing efficiency of the torus. In the models, the outer 
radial extent of the torus scales with the AGN luminosity, and we find the tori to be confined to scales 
less than 5 pc. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - infrared: galaxies 



1. INTRODUCTION 

The unified model fo r active galaxies (|Antonuccill 19931 : 
lUrrv fc Padovanll995[ ) explains the observed differences 
between Type-1 and Type-2 objects as due to orienta- 
tion. Although the unification model may not be uni- 
versally applicable, the most accepted scheme is that 
there is dust surrounding the cen tral region of AGN 
distributed in a toroidal geometry (iPier fc Kroliklll993l: 
Granato fc Danesill994t [Efstathiou fc Rowan- RobinsonI 



19951: iGranato et al.l Il997i r Thus, the central engines 
of Type-1 active galactic nuclei (AGN) can be seen di- 
rectly, resulting in typical spectra with both narrow and 
broad emission lines coming from the Narrow-Line Re- 
gion (NLR) and the Broad-Line Region (BLR), respec- 
tively. In Type-2 AGN, the active nucleus and the BLR 
are obscured by an optically and geometrically thick 
dusty torus, which explains the lack of broad lines in 
their observed spectra. The torus dust grains absorb ul- 
traviolet photons from the central engine and re-radiate 

^ Institute dc Astrofi'sica dc Canarias (lAC), C/Vi'a Lactca, s/n, 
E-38205, La Laguna, Tenerife, Spain; cra@iac.es 

^ Department of Physics and Astronomy, University of Ken- 
tucky, Lexington, KY 40506. 

3 Institute de Estructura de la Materia, CSIC, E-28006 Madrid, 
Spain. 

* Gemini South Observatory, Casilla 603, La Serena, Chile. 

^ Astronomy Department, University of Florida, 211 Bryant 
Space Science Center, P.O. Box 112055, Gainesville, FL 32611- 
2055. 

® Gemini Observatory, Northern Operations Center, 670 North 
Aohoku Place, Hilo, HI 96720. 



them in the IR, peaking at mid-IR wavelengths (7 - 26 

Pioneeri ng work in modellin g dusty tori was per- 
formed by iPier fc KrolikI ()1992f ) . They assumed a uni- 
form dust density distribution of a few parsecs ra- 
dial extent. Ho wever, in order to be tter reproduce 
the observations, I Pier fc KrolikI ()1993D enlarged the 
torus, e.g., to account for the far-IR emission. Thus, 
the coolest material in the torus must be located 
far from the central engine, implying '^100 pc scale 
tori. Further investigations with smooth dust distr i- 
butions were carried out by IGranato fc Danesd ([l99l). 



Granato et al.l 



lEfstiithiou fc Rowan-Robinson' (199E 
(I1997D . and , Siebenm orgcn ct al. (2004a 

However, recent ground-based mid-IR observations of 
nearby Seyferts reveal that the torus s ize is likely re- 
stricte d to a few parsecs. For ex ample, iPackham et al.l 
()2005f) and lRadomski et al.l ()2008D established upper lim- 
its of 2 and 1.6 pc for the outer radii of the Circinus 
galaxy and Centaurus A tori, respectively using data 
from the Gemini Telescopes. Additionally, interferomet- 
ric observations obtained with the MID-infrared Interfer- 
ometric Instrument (MIDI) at the Very Large Telescope 
Interferometer (VLTI) of Circinus, NGC 1068, and Cen- 
taurus A suggest a scenario whe re the torus em ission 
would only extend out to R . = 1 p c (|Tristram et al.ll20G7) , 
R = 1.7 - 2 PC (iJaffe et al.ll2004t|R aban ct al. 2009), and 
R = 0.3 pc (jMeisenheimer et al.ll2007|) . respectively. 

Thus, during the last decade, an intensive search for 
an alternative geometry to the early homogeneous torus 
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models has been c arried out in order to explai n the ob- 
servat ions. Before iPier fc Krolild (|1992l . Il993[ ) publica- 
tions, iKrolik fc BegelmanI (|1988[ ) claimed that smooth 
dust distributions cannot survive within the AGN vicin- 
ity. They proposed that the material in the torus must 
be distributed in a clumpy structure, in order to pre- 
vent the dust grains from being destroyed by the hot 
surrounding gas. The first results of radiative trans- 
fer calculations o f a cl umpy medium were reported 
by iNenkova et all (120021), and further work was done 



bvlDuUemond fc van Bemmell (l2005D.lFritz et al 
lElitzur fc Shlosmanl l 2006( 1 ■ and lBallantvne et al 



Thu s, the c lun ipy dusty to r us mo d els (INenkova et"al 
2001 I2008airbl : iHonig et all l2006t ISchartmann et al 



20081) propose that the dust is distributed in clumps, in- 
stead of homogeneously filling the torus volume. These 
models are making significant progress in accounting for 
the mid-IR emission of AGNs. A fundamental difference 
between clumpy and smooth density distributions of dust 
is that the former implies that both directly-illuminated 
and shadowed cloud faces may exist at different dis- 
tances from the central engine. Thus, the dust temper- 
ature is not a function of the radius only, as was the 
case in the homogeneous models, since illuminated and 
shadowed clouds contribute to the IR emission from all 
viewing angles. Homogeneous torus models had serious 
difficulties in pre dicting the observed differences in the 
variety of SEDs (ISpin oglio fc MalkanI [19891: iFadda et al " 
1998; Kuraszkiewicz e t al.l 120031 : lAl onso-Herrero et al 



2003t iRigbv et all l200l ) and the strengths of the 10 
/xm silicate feature in both Seyfert 1 and Seyfert 2 nu- 



clei (iRoche et al.lll991l: iGranato et"anil997l : iMason et all 
|2006| ). and also the large optical depths along the line 
of sight indicated by X- ray observations (jSimpson et al.l 
119941 : iTozzi et al.l l2006l ). Inhomogeneous torus mod- 
els solve these problems. The clumpy models have 
been employed in several observational studies, e.g.. 
the first Spitzer ana lysis of AGN in the GOODS fields 
(jTreister et al.l |2004|) . and in the analysis of the 10 
/im spatially r esolved spectra of NGC 1068 and NGC 
2110 nuclei (Ma son et a l. 2006, 2009). 

Since the reprocessed radiation from the torus is re- 
emitted in the IR, this range (especially the mid-IR) is 
key to put constraints on the clumpy dusty torus mod- 
els. In comparing the predictions of any torus model 
with observations, the small-scale torus emission must 
be isolated. Large aperture data (i.e., images with an- 
gular resolutions larger than 1-1.5", such as those pro- 
vided by ISO, IRAS, or Spitzer) are strongly contami- 
nated with emission from the host galaxy. For example, 
starburst emission is a.n important source of IR flux in 
the majority of AGN (iNetzer et al. 2007: Ba rmbv et"^ 



20061: iPoUetta et al.ll2007t iRamos Almeida et al.ll2009D . 



Mason et al.l (|2006f ) found for NGC 1068 that all flux 



measurements within apertures smaller than 0.5" were 
well fltted by the models, whereas larger aperture fluxes 
were much higher and presented a different spectral 
shape in the mid-IR. This excess of emission could be due 
to nearby dust o utside the torus, t h at is also emit t ing in 
the IR (see e.g.. lRoche et al.ll2006D . iMason et all ()2006D 
concluded for NGC 1068 that the torus contributes less 
than 30% of the 10 fim flux within apertures larger than 
1". The bulk of this large- aperture flux comes from dust 
in the ionization cone, which occupies a larger volume 



than the torus dust. 

At shorter IR wavelengths (~ 1-2.2 fim) the flux 
contained in apertures of a few arcseconds partly 
comes from the torus and p artly from stellar emission. 
lAlonso-Herrero et al.l (|l996l ) found that for Seyfert 2 
galaxies (Sy2), the fluxes obtained within a 3" aperture 
are dominated by stellar emission in the near-IR, out to ^ 
2 /im. Indeed, iKotilainen et al.l (|1992l ) found for a com- 
plete sample of Seyfert galaxies, most of them Seyfert 1 
(Syl), that starlight contributes significantly to the nu- 
clear emission in the near-IR (J, H, and K bands), even 
within a small projected aperture (of 3"). Consequently, 
high angular resolution observations at these wavelengths 
are crucial to separate the pure nuclear emission from 
that of the host galaxy, which contaminates or hides the 
AGN emission. 

The optimal way to estimate near-IR nuclear fluxes 
uncontaminated by star formation, in as much as possi- 
ble, is the use of either subarcsecond spatial resolution 
data (e.g., from the HST or adaptive optics) or high spa- 
tial resolution ground-based data, extracting the nuclear 
emission in a proper manner (e.g., using PSF subtrac- 
tion). In the mid-IR, it is necessary to isolate the torus 
emission from that produced by dust in the ionization 
cones or dust heated by intense star formation. With 
various high-resolution SEDs covering the near- and mid- 
IR ranges, it is possible to identify SED features that are 
attributable to the torus. 

It is worth mentioning that mid-IR interferometry con- 
stitutes the most precise way of characterizing the torus 
structure using observ ations and compa ring them with 
the torus models (e.g.. iHonig et al]|2006t ). However, the 
signal-to-noise required for VLTI observations strongly 
restricts the number of observable sources to the nearest 
AGN only. 

In this work, we report new high-resolution mid-IR 
imaging data for 18 nearby Seyfert galaxies, for which 
we have estimated nuclear mid-IR fluxes^. We also com- 
piled near-IR high spatial resolution fluxes from the lit- 
erature to construct pure-nuclear SEDs. We have fitted 
th ese AGN SEDs with t he clumpy dusty torus models 
of INenkova et al.l (|2008ai rb[) to constrain the parameters 
that describe the clumpy models when applied to Seyfert 
galaxies. Table [T] summarizes key observational proper- 
ties of the sources in the sample. Section 2 describes the 
observations, data reduction, and compilation of near-IR 
fluxes. Section 3 presents the main observational results, 
and §4 presents the modelling results. We draw conclu- 
sions about the clumpy torus models and AGN obscura- 
tion in general in §5. Finally, Section 6 summarizes the 
main conclusions of this work. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Mid-IR Imaging Observations 

Ground-based mid-IR high-angular resolution obser- 
vations of 18 nearby active galaxies were carried out 
over the past years for a variety of science drivers. We 
make use of this archive of data in this paper. Most of 
these sources are Type-2 Seyferts, but the sample also 

^ The mid-IR images of Circin us, NGC 4151 , and Cen- 
taurus A are a l ready publi shed by IPackham et al.l I I2005I) and 



IRadomski et all l|2003l . 120081 '). respectively. However, we present 
new calculations of the nuclear mid-IR fluxes in this work. 
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TABLE 1 
Basic Galaxy Data 



n 1 








i_^lk5 Lfllice 


Scalp 










(^ivipc; 


f pc arcscc ) 


Centaurus A 


Sy2 


so pec 


0.0018 


3.5 ^ 


17 


Circinus 


Sy2 


SA(s)b; 


0.0014 


4 a 


20 


IC 5063 


Sy2 


SA{s)0+: 


0.0113 


45 


219 


Mrk 573 


Sy2 b 


{R)SAB{rs)0+: 


0.0172 


69 


334 


NGC 1386 


Sy2 


SB{s)0+ 


0.0029 


11 


56 


NGC 1808 


Sy2 


{R'l)SAB{s:)b 


0.0033 


11 =1 


64 


NGC 3081 


Sy2 


{Ri)SAB(r)0/a 


0.0080 


32 


155 


NGC 3281 


Sy2 


SAB(rs+)a 


0.0107 


43 


208 


NGC 4388 


Sy2 


SA{s)b: sp 


0.0084 


34 


163 


NGC 5728 


HII/Sy2 


(Ri)SAB(r)a 


0.0094 


38 


182 


NGC 7172 


HII/Sy2 


Sa pec sp 


0.0087 


35 


169 


NGC 7582 


Sy2 


(R'i)SB{s)ab 


0.0053 


21 


103 


NGC 1365 


HII/Syl.8 


(R')SBb(s)b 


0.0055 


18 


107 


NGC 2992 


Syl.9 


Sa pec 


0.0077 


31 


149 


NGC 5506 


Syl.9 ^ 


Sa pec sp 


0.0062 


25 


120 


NGC 3227 


Syl.5 


SAB(s) pec 


0.0039 


17 


82 


NGC 4151 


Syl.5 


(R')SAB(rs)ab: 


0.0033 


13 


64 


NGC 1566 


Syl 


(R'i)SAB(rs) 


0.0050 


20 


97 



Note. — Classification, morphological class, and spectroscopic redshift of our targets taken from the NASA/IPAC Extragalactic 
Database (NED). The distance to the sources and the physical scale have been obtained using Ho=75 km s~^ Mpc~^. 

For the cases of the most nearby objects Centaurus A, Circinus, NGC 1386, NGC 180 8 , NGC 1365, NGC 3 2 27, and NGC 
4151 the values of the distance to the galaxies have been taken from [Radomski et al.I (12008'), 'Packham et al.l II2005I ). IBennert et al.l 
(|200i), IJimenez-Bailon et al.l lHooD, ISilbcrma nn et al.l 119991) . IDavies et a l. (2006), and Radomski ct al. (2 003), respectively. ^ Mrk 573 
and NGC 5506 have rece ntly been classified as obscured Narrow-line Syl galaxies ( iNagar et al.i i2002l ; IDewanean &: Griffithsl 120051 : 
IRamos Almeida et ani2008l) . 



includes two Seyfert 1.9 (Syl.9), one Seyfert 1.8 (Syl. 8), 
two Seyfert 1.5 (Syl.5). and one Syl galaxy (see Table 

The first set of observations was obtained with the 
University of Florida mid-IR camera/spectrometer OS- 
CIR, in 1998 December at the 4 m Blanco Telescope at 
Cerro Tololo Inter-American Observatory (CTIO) and in 
2001 May at the Gemini North Telescope. This instru- 
ment employs a Rockwell 128x128 pixel Si:As blocked 
impunity band (BIB) detector array, optimized for the 
wavelength range between 8 and 25 /Ltm. On CTIO, 
the pixel scale is 0.183", with a field of view (FOV) of 
23.4" X 23.4", and the resolutions obtained were ^1" at 
both 10.7 and 18.2 /xm. On Gemini North, OSCIR has 
a plate scale of 0.089" pixel" ^ and a total FOV of 11" x 
11". The achieved resolutions were ^^0.5" at 10.7 /im and 
~0.6" at 18.2 /im. Another set of observations was per- 
formed with the mid-IR camera/spcctrog raph T-ReCS 
(Th ermal- Region Camera Spectrograph; iTelesco et al.l 
[l998i) on the Gemini-South telescope. T-ReCS uses a 
Raytheon 320x240 pixel Si:As IBC array, providing a 
plate scale of 0.089" pixel" ^, corresponding to a FOV 
of 28.5" X 21.4". The resolutions obtained were 0.3- 
0.4" at 8.8 /im, -0.5" at 10.4 /tm, and 0.5-0.6" at 18.3 
/im. The last set of observations w as made with the m id- 
IR camera/spectrograph Michelle (jClasse et al.lll997l ) on 
the Gemini North Telescope. The Michelle detector is a 
Si:As IBC array with a format of 320x240 pixels. Con- 
figured in its imaging mode, Michelle has a 0.10" pixel"^ 
plate scale, and the resolutions of our observations were 
~0.4" at 11.3 /im and -0.5" at 18.1 /im. This plate 
scale translates to a FOV of 32" x 24". A summary of 
the observations is reported in Tabled 

The standard chopping-nodding technique was used to 
remove the time-variable sky background, the telescope 
thermal emission, and the so-called 1/f detector noise. 



The chopping throw was 15", and the telescope was nod- 
ded every 30 s. All data were reduced using in-house- 
developed IDL routines. The difference for each chopped 
pair for each given nodding set was calculated, and the 
nod sets were then differenced and combined until a sin- 
gle image was created. Chopped pairs obviously compro- 
mised by cirrus, high electronic noise, or other problems 
were excluded from any further reduction process. 

Observations of flux standard stars were made for the 
flux calibration of each galaxy through the same filters. 
The uncertainties in the flux calibration were found to 
be -5-10% at N band and -15-20% at Qa band. PSF 
star observations were also made immediately prior to or 
after each galaxy observation to accurately sample the 
image quality. These images were employed to determine 
the unresolved (i.e., nuclear) component of each galaxy. 
The PSF star, scaled to the peak of the galaxy emission, 
represents the maximum contribution of the unresolved 
source, where we integrate flux within an aperture of 
2" for the Gemini data, and of 4" for the CTIO data 
(NGC 1365, NGC 1808, and NGC 7582). The residual 
of the total emission minus the scaled PSF represents 
the host galaxy contribution. A detailed study of the 
extended near nuclear structures of the galaxies in our 
sample will be the subject of a forthcoming paper (J. 
T. Radomski et al. 2009, in preparation). We require a 
flat profile in the residual (for a realistic galaxy profile) 
and therefore reduce the scale of the PSF from match- 
ing the peak of the galaxy emission (100%) to obtain the 
unresolved fluxes reported in Table [3] Figure [1] shows 
an example of PSF subtraction at various levels (in con- 
tours) for the Sy2 galaxy NGC 4388 in the N' Mi chelle 
filter, perfo r med f ollowing I Radomski et al.l ()2002f ) and 
ISoifer et all (|2003f ). The residual profiles fr om the dif- 
ferent scales demonstrate the best-fitting result, scaled 
to 90% of the peak flux. The uncertainty in the un- 
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TABLE 2 
Summary of Mid-IR Observations 



Galaxy 


Filters 


Instrument 


Telescope 


Observation 


On-Source Time (s) 


PSF FWHM 










epoch 


N band 


Q band 


N band 


Q band 


Ccntaurus A 


Si2, Qa 


T-ReCS 


Gemini S 


Jan 2004 


2000 


1550 


0.30 


0.53 


Circinus 


Si2 Oa 


T-ReCS 


Gemini S 


Feb 2004 


109 


109 


0.33" 


U.DO 


IC 5063 


Si2, Qa 


T-ReCS 


Gemini S 


Jul 2005 


1 "^n 

-LOU 


ou^ 






Mrk 573 


N, Qa 


T-ReCS 


Gemini S 


Dec 2003 


01 T 
Zl ( 


Zi ( 


U.oD 


U.04 


NGC 1386 


N, Qa 


T-ReCS 


Gemini S 


Dec 2003 


217 


217 


0.31" 


0.54" 


NGC 1808 


N, IHW18 


OSCIR 


CTIO 4m 


Dec 1998 


300 


300 


0.94" 


1.02" 


NGC 3081 


Si2, Qa 


T-ReCS 


Gemini S 


Jan 2006 


130 


304 


0.30" 


0.56" 


NGC 3281 


N, Qa 


T-ReCS 


Gemini S 


Jan 2004 


260 


455 


0.34" 


0.58" 


NGC 4388 


N', Qa 


Michelle 


Gemini N 


May 2006 


549 


733 


0.34" 


0.50" 


NGC 5728 


Si2, Qa 


T-ReCS 


Gemini S 


Jul 2005 


130 


304 


0.35" 


0.56" 


NGC 7172 


N 


T-ReCS 


Gemini S 


May 2004 


305 




0.51" 




NGC 7582 


N, IHW18 


OSCIR 


CTIO 4m 


Dec 1998 


250 


250 


0.76" 


0.99" 


NGC 1365 


N, IHW18 


OSCIR 


CTIO 4m 


Dec 1998 


482 


482 


0.92" 


1.03" 


NGC 2992 


N', Qa 


Michelle 


Gemini N 


May 2006 


730 


1095 


0.32" 


0.53" 


NGC 5506 


N', Qa 


Michelle 


Gemini N 


Apr 2006 


546 


729 


0.36" 


0.51" 


NGC 3227 


N' 


Michelle 


Gemini N 


Apr 2006 


300 




0.39" 




NGC 4151 


N, IHW18 


OSCIR 


Gemini N 


May 2001 


360 


480 


0.53" 


0.58" 


NGC 1566 


Si2, Qa 


T-ReCS 


Gemini S 


Sep 2005 


152 


304 


0.30" 


0.53" 



Note. — Images were obtained in the 8.74 fim (Si2, AA = 0.78 fira at 50% cut-on/off), 10.36 /^m (N, AA = 5.27 fim), and 18.33 fim (Qa, 
AA = 1.5 fj.m) filters with T-ReCS; in the 11.29 ^tm (N', AA = 2.4 fim) and 18.11 ^tm (Qa, AA = 1.9 ^im) filters with Michelle; and in the 
10.75 fim (N, AA = 5.2 ^tm) and 18.17 ^tm {IHW18, AA = 1.7 fim) filters with OSCIR. 

resolved fluxes determination from PSF subtraction is 
~10-15%. Thus, we estimated the errors in the flux den- 
sities reported in Tablc[3]by adding quadratically the flux 
calibration and PSF subtraction uncertainties, resulting 
in errors of ~15% at N band and of ^-^25% at Qa band. 

For the galaxies that clearly appear spatially unre- 
solved in our mid-IR images (Mrk 573, NGC 1566, and 
NGC 5728), aperture photometry was employed to de- 
rive their nuclear density fluxes (reported in Table [3]) . 
The rule for their calculation was to ensure that all the 
unresolved emission was collected in the chosen aperture, 
with the minimum contribution of foreground contribu- 
tion as possible. Aperture correction was applied to these 
fluxes, employing the images of the corresponding stan- 
dard stars in both N and Qa filters. 



2.2. Compilation of Near-IR High Spatial Resolution 

Data 

To construct well-sampled IR SEDs, we compiled near- 
IR high spatial resolution nuclear fluxes from the liter- 
ature (reported in Table d]). These fluxes correspond 
to the observed emission from the nuclear region of the 
galaxies (unresolved component), with typical spatial 
scales of the same order as those of our nuclear mid-IR 
fluxes. NICMOS/HST fluxes are available for 10 galax- 
ies, having FWHM for an unresolved PSF ~0.13" in the 
F160W flher and ~0.09" in FllOW. For the cases of Ccn- 
taurus A and Circinus, there are also diffraction-limited 
and near-diffraction- limited adaptive optics NACO /VLT 
fluxes, respectively. (The FWHM of Centaurus A un- 
resolved component is 0.1" in the near-IR bands and 
for Circinus 0.2" in the near-IR bands and < 0.13" in 
the L' and M' bands.) We also compiled seeing-limited 
near-IR nuclear fluxes from the literature, obtained with 
data from IRCAM3 on the 3.8 m United Kingdom IR 
Telescope (UKIRT), with NSFCam at the 3 m NASA 
IRTF telescope, with the IRAC-1 IR array camera on 
the ESQ 2.2 m telescope, and with ISAAC at the 8 
m VLT. The angular resolution of these natural seeing- 



TABLE 3 
Unresolved Mid-IR Fluxes 



Galaxy Level of PSF subtraction Flux Density (mjy) 





N band 


Qa band 


N band 


Qa band 


Centaurus A 100% 


100% 


710 


2630 


Circinus 


90% 


90% 


5620 


12790 


IC 5063 


100% 


100% 


399 


2790 


Mrk 573=* 


100% 


100% 


177 


415 


NGC 1386 


70% 


70% 


147 


457 


NGC 1808 


40% 


60% 


227 


560 


NGC 3081 


100% 


50% 


83 


231 


NGC 3281 


80% 


100% 


355 


1110 


NGC 4388 


90% 


100% 


195 


803 


NGC 5728=^ 


100% 


100% 


25 


184 


NGC 7172 


90% 




68 




NGC 7582 


70% 


50% 


195 


527 


NGC 1365 


90% 


100% 


321 


640 


NGC 2992 


90% 


90% 


175 


521 


NGC 5506 


100% 


100% 


873 


2200 


NGC 3227 


100% 




401 




NGC 4151 


90% 


100% 


1320 


3200 


NGC 1566" 


100% 


100% 


29 


117 


Note. — 


The percentaj 


;es of PSF subtraction level 


are reported 



in the employed filters (listed in Table [2J. Errors in fiux densities 
are dominated in general by uncertainties in the flux calibration 
and PSF subtraction (~15% at N band and ~25% at Qa band). 

For the galaxies Mrk 573, NGC 1566, and NGC 5728, aper- 
ture photometry was performed with the following aperture radii: 
0.89", 0.45", and 0.62", respectively. These fluxes include the cor- 
responding aperture corrections. 

limited images is 0.6-0.7" in K band and 0.6-0.9" in L 
band. These ground-based near-IR nuclear fluxes were 
estimated either using aperture photometry (scaling an 
annulus of the J-band images to the observed co unts in H 
and K, ISimpsonlfl998llGalhano fc Alloinll2008D. surface 
brightness pr oflle deconvolution ( Quillen c t al. 200ll), or 
PSF scaling (|Alonso-Herrero et al.ii2001i. i2003il . Due to 
the lack of high spatial resolution nuclear fluxes of Seyfert 
galaxies in the literature, the SEDs in our sample range 
from very well-sampled SEDs (e.g., Circinus and Mrk 
573) to galaxies for which there are no high spatial res- 
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olution ncar-IR nuclear fluxes available (NGC 1566 and 
NGC 5728). 



3.1. 



Fig. 1. — 11.29 fim contour plots of NGC 4388, the PSF star, and the scaled subtraction of the PSF for this galaxy at various levels 
(100%, 90%, and 80%). The residuals of the subtraction in the lower right panel show the host galaxy profile, with the residual at 90% 
scale representing the best fit. 

spatial resolution data, we would find umir = oliras- 
However, that is not the case, and except for a couple 
of galaxies for which both slopes are similar (Mrk 573, 
NGC 1365, and NGC 1808), the rest of sources are far 
from the umir — oliras relationship. 

An analogous comparison between our unresolved mid- 
IR fluxes and the IRAS measurements comes from the 
N/IRAS 12 /im and Q/IRAS 25 ^lm. ratios (see Tabic [5]). 
The low values of these ratios indicate the large amount 
of emission related to stellar processes that is contribut- 
ing to the IRAS fluxes in contrast to our nuclear mid-IR 
fluxes, that are representative of torus emission. The 
lowest values of the ratios correspond to Centaurus A, 
NGC 1365, NGC 1808, NGC 1566, and NGC 7582. 

In summary, the high spatial resolution mid-IR mea- 
surements provide a spectral shape of the SEDs that is 
different from that of large aperture data SEDs. The nu- 
clear SEDs presented in this work will allow to us charac- 
terize the torus emission and consequently, to use torus 
models to constrain the distribution of dust in the im- 
mediate AGN vicinity. 



3. SED OBSERVATIONAL PROPERTIES 

High versus Low Spatial Resolution Spectral Energy 
Distributions 

Important observational constraints for torus mod- 
elling arise from the shape of the IR SEDs of Seyfert 
galaxies, as the bulk of the torus emission is concentrated 
in the IR range, peaking at mid-IR wavelengths (7 - 26 
/zm). When large aperture data (i.e., of angular resolu- 
tions larger than 1-1.5") such as those from ISO, Spitzer, 
or IRAS, are employed for constructing SEDs, the IR 
fluxes are a mixture of AGN plus host galaxy emission. 
In particular, for the case of Sy2 galaxies, stellar emis- 
sion dominates the n ear-IR fluxes at A 5, 2.2 /i m even 
within 3" apertures (jAlonso-Herrero et al. I [19961 ). Even 
for Syl galaxies, where the direct emission of the AGN is 
generally the dominant contribution to the near-IR flux, 
iKotilainen et al.l (|l992l ) found that starlight contributes 
significantly to their nuclear emission. The stellar con- 
tribution can also be significant in large-aperture mid- 
IR me asurements. For example, iRamos Almeida et al.l 
(|2007t) found that the ISOCAM mid-IR fluxes (with an 
effective resolution in the diffraction limit of ^4") of a 
sample of ^60 Seyfert galaxies are dominated by AGN 
emission in Syl, although there is a non- negligible con- 
tribution from the host galaxy that becomes larger in 
Sy2. 

Using the data described in Section [2] we have con- 
structed high spatial resolution SEDs from ~ 1-18 /im. 
In order to compare them with large aperture measure- 
ments, we have compiled IRAS 12 and 25 /^m fluxes for 
the galaxies in our sample, taken from NED (Table [5|). 
We compare the nuclear mid-IR spectral slopes calcu- 
lated from our N and Q data points {(Xmir] column 
3) with the slopes from the IRAS 12 and 25 /im mea- 
surements [aiRAS] column 4). If the spectral shape of 
the IRAS SEDs were similar to that of our mid-IR high- 



3.2. Average Seyfert 2 Spectral Energy Distribution 

In order to derive general properties of Sy2 galaxies, 
which constitute the majority of our sample, we con- 
structed an average Sy2 SED, considering only the high- 
est angular resolution data to avoid as much as possi- 
ble the stellar contamination. This average Sy2 tem- 
plate will be used to compare with the individual SEDs 
of the Sy2 analyzed in this work, and also with the 
intermediate-type Seyferts. Out of the 18 galaxies in to- 
tal (see Table[T|), there are 10 "pure" Type-2 Seyferts and 
two HII/Sy2 composite objects. The rest of the sample 
consists of other Seyfert types (Syl. 9, Syl. 8, Syl. 5, and 
Syl). For constructing the Sy2 average SED, we first ex- 
cluded the composite objects and the intermediate-type 
Seyferts. Second, to ensure that we are collecting the 
best-measured non-stellar nuclear emission, we exclude 
galaxies with mid-IR data from the 4 m CTIO, because 
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TABLE 4 

High Spatial Resolution Near-IR Fluxes 



Galaxy 




Flux Density (: 


mJy) 




Filters 


Reference(s) 




J band 


H band 


K band 


L band 


M band 






v_>6ntaiiriis A 


i .omu. i 


4.5±0.3 


34±2 


zuu^^u 






a 


Circinus^ 


1.6±0.2 


4.8±0.7 


19±2 


o o n 1 o o 


1900±190 


1160W, NACO J,K,L,M 


b 


bubo 




0.3±0.1 








b 160W 


c 


A/Trlf ^7'^ 

iV±l rv ij 1 O 


KJ . -LO^U. UU 


0.54±0.04 


3.2±0.6 


1 Q Q-L'i Q 
±0.0^.0. o 


41.3i8.3 


Fl 1 OW Fl 60W NSFCam K' T, M 


(J 


INCjC lOOD 




0.2±0.1 








t ioUW 


c 




15.5±4.5 




30.5±8.5 






TSA AC T Ks T,' 




NGC 3081 




0.22±0.13 








F160W 


c 


NGC 3281 




1.3±0.2 


7.7±0.8 


103±9 


207±25 


IRAC-1 H,K, IRCAM3 L',M 


f 


NGC 4388 


0.06±0.02 


0.71±0.28 




40±8 




F110W,F160W, NSFCam L 


d 


NGC 5728 
















NGC 7172 




<0.4 


3.4±0.7 


30±6 


61±12 


IRCAM3 H,K,L',M 


g 


NGC 7582 




22.6±2.3 




117±18 


142±21 


F160W, ISAAC L,M 


c,h 


NGC 1365 




8.3±0.8 


78±8 


205±41 


177±35 


F160W 


1 


NGC 2992 




<1 


2.8±0.6 


22.7±4.5 


35.7±7.1 


IRCAM3 H,K,L',M 


g 


NGC 5506 


13.8±2.8 


59±12 


120±24 


340±68 


530±106 


IRCAM3 J,H,K,L',M 


g 


NGC 3227 




10.6±1.1 


22.6±4.5 


46.7±9.3 




P160W, NSFCam K',L 


d 


NGC 4151 


69±7 


104±10 


177±35 


325±65 


449±34 


F110W,F160W, NSFCam K',L, IRCAM3 M 


d 


NGC 1566 



References. - — (a ) iMciscnlicimc r ct aLI r2007f): ( b) Pricto et al. (2004); (c) Quillcn ct al. ( 20011 ) ; (d) Alonso-Hcrrcro ct al.l II2003D : (e) 
IGalliano fc AlloinI (I2008D : (f)lSimpsoii ( 1998D ; (g) lAlonso-H crrcro ct al. (2001): (li) Pricto ct al. (2003); (i) ICarollo ct a l. (2002) 

Note. — Ground-based instruments and telescopes are: NACO and ISAAC on the 8 m VLT, NSFCam on the 3 m NASA IRTF, IRCAM3 on 
the 3.8 m UKIRT, and IRAC-1 on the 2.2 m ESQ telescope. Measureme nts in the Fl l OW a nd F160W filters are from NICMOS on HST. 
^ In the case of Circinus, we also use the 2.42 fim flux (31±3 mJy) from lPrieto et al.l II2004D . 



TABLE 5 

Spectral Shape Information 



Galaxy 








OtIRAS 


H/N 


N/Q 


N/IRAS 12 


Q/IRAS 25 


Average 


3.1±0.9 


3.6±0.8 


2.0±0.2 




0.003 


0.23 






Centaurus A 


2.8±0.8 


3.4±0.8 


1.8±0.2 


0.3±0.1 


0.006 


0.27 


0.03 


0.09 


Circinus 


3.6±1.0 


4.5±1.0 


l.litO.l 


1.8±0.2 


0.001 


0.44 


0.30 


0.19 


IC5063 


3.8±1.5 


4.2±1.7 


2.7±0.4 


1.8±0.2 


0.001 


0.14 


0.36 


0.65 


Mrk 573 


2.9±0.8 


3.4±0.8 


1.7±0.1 


1.5±0.1 


0.003 


0.43 


0.63 


0.49 


NGC 1386 


3.3±1.2 


3.7±1.4 


2.0±0.2 


1.4±0.1 


0.001 


0.32 


0.28 


0.31 


NGC 1808 


1.8±0.3 


1.7±0.2 


1.7±0.1 


1.6±0.1 




0.41 


0.04 


0.03 


NGC 3081 


3.0±1.1 


3.5±1.3 


1.4±0.1 




0.003 


0.36 






NGC 3281 


2.6±0.6 


3.1±0.7 


2.0±0.2 


1.4±0.1 


0.004 


0.32 


0.39 


0.42 


NGC 4388 


3.0±0.8 


3.3±0.7 


2.7±0.4 


1.7±0.2 


0.004 


0.24 


0.19 


0.22 


NGC 5728 






2.7±0.2 


2.2±0.3 




0.14 


0.12 


0.18 


NGC 7172 


1.8±0.2 


1.8±0.2 




l.OiO.l 


<0.006 




0.16 




NGC 7582 


1.2±0.1 


l.litO.l 


1.9±0.1 


1.6±0.1 


0.116 


0.37 


0.08 


0.07 


NGC 1365 


1.8±0.4 


2.0±0.5 


1.3±0.1 


1.4±0.1 


0.026 


0.50 


0.06 


0.04 


NGC 2992 


2.4±0.5 


2.7±0.4 


2.1±0.2 


l.OiO.l 


<0.006 


0.34 


0.24 


0.33 


NGC 5506 


1.7±0.3 


2.0±0.3 


1.2±0.1 


1.6±0.1 


0.068 


0.40 


0.68 


0.53 


NGC 3227 


1.8±0.2 


1.8±0.2 




0.9±0.1 


0.026 




0.43 




NGC 4151 


1.4±0.2 


1.3±0.1 


1.6±0.1 


1.2±0.1 


0.079 


0.41 


0.66 


0.66 


NGC 1566 






2.0±0.2 


0.6±0.1 




0.25 


0.02 


0.04 



Note. — Fitted spectral indexes (fi^ a in the whole range {aiR, from ~1 to 18 /im), in the near-IR (ajv/n, from ~1 to ~9 fj,m), 

and in the mid-IR using the N and Q band data points). The IRAS slopes {ctmASy determined using the 12 and 25 fim IRAS 

fluxes) are shown also for comparison. H/N and N/Q band ratios are reported, together with the N/12 ^m and Q/25 ^m IRAS ratios. 

the average template: Centaurus A, Circinus, IC 5063, 
Mrk 5738, NGC 1386, NGC 3081, and NGC 4388. Their 
IR SEDs are shown in Figured We normalize the SEDs 
at 8.8 /im in order to emphasize the comparison between 
the N- and both the H- and Q-band data points, that 
provides useful information about the obscuration (see 
Section ISTT]) . The measured N/Q ratio for the average 
template is 0.23±0.14, which is representative of the Sy2 
in our sample. 

^ Note that Mrk 573 has been recently reclassifled by 
IRamos Almeida et al.l (|2008l ) as an obscured NLSyl, but the IR 
SED remains similar to that of the ordinary Sy2 galaxies. 



of their lower spatial resolution compared with that of 
the 8 m Gemini telescopes. This restriction excludes 
NGC 1808 and NGC 7582. Indeed, these two galaxies 
are undergoing intense episodes of nuclear star forma- 
tion, making it more difficult to isolate the AGN emission 
from that of the surrounding host galaxy. Finally, in the 
construction of the average template, we use only HST 
(NICMOS) and VLT (NACO) adaptive optics fluxes in 
the near-IR. This excludes NGC 3281, for which all the 
compiled data are ground-based, and a few points in the 
SEDs of some of the other galaxies (e.g., Mrk 573 and 
NGC 4388). 

In the end, we have seven pure Sy2 for constructing 
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Thus, the average template is contructcd with near- 
and mid-IR data having resolutions < 0.6", to ensure the 
isolation of the torus emission. The best-sampled SED, 
Circinus, defines the wavelength grid, and we interpo- 
lated nearby measurements of the other six galaxies onto 
this scale. We did not interpolate the sparse observations 
of IC 5063 and NGC 3081. The interpolated fluxes were 
used solely for the purpose of deriving the average Sy2 
SED template, shown in Figure [21 The error bars corre- 
spond to the standard deviation of each averaged point, 
except for the 8.8 /xm point (the wavelength chosen for 
the normalization). In this case, we assigned a 15% er- 
ror, the nominal percentage considered for the N-band 
flux measurements. 

3.3. Type-2 Seyfert SEDs 

The shape of the Sy2 SEDs in our sample is gener- 
ally very steep, compared with the intermediate-type 
Seyferts. In geiieral, S y2 have steeper 1-10 /xm SEDs 
than S vl (iRiea [TOTSl: fedclson et al.l[T987l: iWard et all 
1987t iFadda et al.nT998; Alonso-Herrero et all 120011 



20031) . Nevertheless, we fin d a variety of spectral shape s 



for Sy2, in agreement with lAlonso-Herrero et al.l (|2003f ). 

We measured the 1-18 iim IR slopes (fj^ a of the 

individual Sy2 galaxies and of the average template (Ta- 
blets]). The IR slopes are representative of the whole SED 
shape. We also calculated the near-IR (ajv//?, from '^l to 
~9 fim) and mid-IR spectral indexes {aMiR-, using the N 
and Q band data points) to compare them. A flat near- 
IR slope indicates an important contribution of the hot 
dust emission (up to ^^1000-1200 K: iRickc &: Lebofskvi 
119811: lBarvainij[l987l) that comes from the inner part of 
the torus, and consequently, must be hidden in smooth 
torus descriptions of Sy2. 

The Sy2 SEDs are typically steep through the IR, with 
<^IR = 3.1 ± 0.9 for the average template. However, 
the individual galaxies do show some variety, with a/^ 
ranging from 1.8 to 3.8, being even shallower in the un- 
certain NGC 7582. Our re sults are similar to those of 
lAlonso-Herrero et al.l (|2003f ). who found afj^"''° > 2.6 
for the 1 to 16 fim slopes of the Sy2 galaxies in the 
expanded CfA sample. Comparing the near- and mid- 
IR slopes for the Sy2 reported in Table [5l we find 
that ajv/ij > ctMiR for all of them, except for NGC 
7582. For the average template, a^iR = 3.6 ± 0.8 and 
aMiR = 2.0 ± 0.2. 

The flatter slopes arise naturally in some configura- 
tions of clumpy torus models with some near-IR contri- 
bution from the hot dust faces of directly-illuminated 
clouds, even when the central engine is blocked from 
view. In contrast, smooth torus models do not gener- 
ally produce a range of IR spectral slopes. Indeed, large 
opacity torus mo dels without extended cone component 
predi ct am > 4 (jPier fc Krohklll993tlGranato fc Danesd 
[l99l . Homo gcneus t ori with lower optical depth can 
flatte n the IR slopes (jGranato et al.l Il997l : iFadda et al.l 
|1998[) . but they predict the 10 /xm sihcate feature in 
strong emission for Typ e-1 nuclei, which is not ob- 
served (jRoche et al.lll991t). The extended conical com- 
ponen t included in the Efstathiou fc Rowan-RobinsonI 
torus models flattens the 3-5 /xm continuum of 
Sy2. However, the high resolution measurements consid- 
ered here should be mostly free of cone contamination. 
Thus, clumpy dusty torus models seem to be the most 



promising way for reproducing the observed range of Sy2 
SEDs. 



3.4. Intermediate- Type Seyfert SEDs 

The Syl. 5 , Syl .8, and Syl.9 galaxies were added by 
lOsterbrockl ()1981[ ) to the Seyfert classiflcation scheme 
to account for those galaxies with weaker featureless 
continua, larger broad Balmer decrements, and weaker 
broad components of the p ermitted lines than Type-1 
Seyferts. [OstcrbrocS (|l98lh suggested that the observed 
characteristics of these intermediate Seyferts could be 
due to the red dening of both the continuu m and the 
BLR by dust. lAlonso-Herrero et all (|2003D found for 
the CfA Seyfert sample, that the majority of galax- 
ies optically classified as Syl.8 and Syl.9 display cither 
IR spectral indexes and SEDs similar to those of Syl 
or intermediate between Sy l and Sy2. For the Syl.5, 
lAlonso-Herrero et al.l (|2003[) found the SEDs to be prac- 
tically identical to those of Syl. 

Figure O shows the SEDs of the five intermediate-type 
Seyferts considered in this work (two Syl.5, one Syl.8, 
and two Syl.9), together with the mean Sy2 template for 
comparison, interpolated onto the mean Sy2 wavelength 
grid. Spectral indexes are reported in Tabled The spec- 
tral slopes of the intermediate- type Seyferts are shallower 
than those of the Sy2 (mean slope ajji ~ 2.0 ± 0.4 for 
Syl.8 and Syl.9 an d ajR = 1.6 ± 0.3 for Syl.5 , as op- 
posed to 3.1 ±0.9). \M onso-Herrero et ah! ()2003D report 
IR spectral indices (from 1 to 16 /xm) of types 1.8 and 
1.9 similar to our measurements {af^"^° = 1.8-2.6), and 
within the interval a^^°""° = 1.5-1.6 for Syl and Syl.5 
galaxies. The steepest SED among the intermediate-type 
Seyferts is in the Syl.9 NGC 2992, which is more similar 
in shape to the Sy2 SEDs. This result could be due to 
an incorr ect or genuinely cha nged classification of this 
galaxy, as iTrippe et al.l (|2008f ) claim (see Appendix E| • 

The near-IR slopes of the Syl.8 and Syl.9 have in- 
termediate values between those of Sy2 and Syl.5, with 
(Xnir. > ctMiR- For the Syl.5 galaxy NGC 4151, aAT/i? < 
oiMiR, indicating the importance of the hot dust contri- 
bution in the SEDs of this type of nucleus. 

The H/N and N/Q fiux ratios show comparable differ- 
ences among Seyfert types (Table [5]). The H/N ratio is 
larger for Syl.5 (0.05±0.04) than for Sy2 (0.003±0.002 
for the average template). The N/Q ratios of the 
intermediate-type Seyferts are also larger than those of 
Sy2 (mean values of 0.41±0.08 for Syl.8 and Syl.9 and 
0.41±0.06 for Syl.5, as opposed to 0.23±0.14). 

In summary, the slope of the IR SED is in general cor- 
related with the Seyfert type. Sy2 show steeper SEDs, 
and intermediate-type Seyferts are flatter. Seyferts 1.8 
and 1.9 present intermediate values of the IR slope and 
H/N ratio between Sy2 and Syl.5. However, we find 
a range of spectral shapes among the Sy2 galaxies, and 
some intermediate-type SEDs have the same slopes as the 
Sy2. This cannot be reconciled with the predictions of 
early torus models, since large optical depth homogenous 
tori strictly predict steep SEDs for Sy2 and fiat SEDs 
for Syl. We do not observe such a strong dichotomy 
(see also lAlonso-Herrero et al]l2003| ). Torus models us- 
ing lower opacities are observationally discarded ( §3.3p 
and models including an extended conical component do 
not seem adecuate for being applied to our high spatial 
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Fig. 3. — Observed IR SEDs normalized at 8.8 fim for the Syl.5, Syl.8, and Syl.9 galaxies in our sample (in color and different symbols). 
The average Sy2 SED is displayed for comparison, displaced in the Y-axis for clarity. Note the difference in steepness between the Sy2 
mean SED and the intermediate types, with the exception of the Syl.9 NGC 2992. 



resolution data . Thus , we pursue clumpy torus models of 
iNenkova et al.l (120021 ) to reproduce the observed SEDs. 

4. SED MODELLING 



4.1. 



Clumpy Dusty Torus Models 

The clumpy dusty torus models of INenkova et al.l 
(|2002, 2008a,b) propose that the dust surrounding the 
central engine of an AGN is distributed in clumps, in- 
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stead of homogeneously filling the torus volume. These 
clumps are distributed with a radial extent characterized 
by the parameter Y = Ro/Rd, where Ro and Rd are the 
outer and inner radius of the toroidal distribution, re- 
spectively (see Figure |4]). The inner radius is defined by 
the dust sublimation temperature {Taub ~ 1500 K), with 
Rd = 0.4 (1500 K r^-ij2.6(^/j^Q4.5 j,j.g s-i)0-5 pc. Within 
this geometry, each clump has the same optical depth 
{tv, defined at V). The average number of clouds along 
a radial equatorial ray is A'o . The radial density profile is 
a power-law (oc r~''). A width parameter, a, character- 
izes the angular distribution, which has a smooth edge. 
The number of clouds along the line of sight (LOS) at an 
inclination angle i is A^los(*) = -^o e^"^*"^"' \ 



LOS 
■4 



o o " o <^ 



O o <^ <^ 



O Q O O o 



R 1 ^ ^ 



.0 



Fig. 4 . — Scheme of the clumpy torus described in lNenkova et all 
l|2008al lH). The radial extent of the torus is defined by the outer 
radius (Ro) and the dust sublimation radius (Rd)- All the clouds 
are supposed to have the same ty, and a characterizes the width 
of the angular distribution. The number of cloud encounters is 
function of the viewing angle, i. 

The radiative transfer equat ions are so lved for each 
clump using the DUSTY code (jlvezic et a l. 1999), with 
solutions depending mainly on the location of each clump 
within the torus, its optical depth, and the chosen dust 
composition. We adopt a dust extinction profile cor- 
responding to the OH Mc dust (i.e., the standard cold 
oxygen-rich ISM dust of lOssenkopf et al.|[r992:) . DUSTY 
includes dust absorption, emission, and scattering com- 
ponents. The total torus emission is calculated by in- 
tegrating the source function of the whole number of 
clumps convolved with the radiation propag ation proba- 
bility along the torus (jNenkova et al.l [20021 ). The direct 
AGN emission may also be included in the resulting SED, 
which is appropriate for Type-1 nuclei and some inter- 
mediate types. 

4.2. BayesClumpy 

The clumpy database contains ~10^ models, calcu- 
lated for a fine grid of model parameters. The inherent 
degeneracy between model parameters has to be taken 
into account when fitting the obs ervables. To this aim, 
lAsensio Ramos fc Ramos Almeid a (2009) recently devel- 
oped a Bayesian inference tool (BayesClumpy) that 
extracts as much information as possible from the ob- 
servations. They applied interpolation methods that are 



able to derive models for different combinations of the six 
clumpy model parameters described above even if they 
are not present in the original database. 

The synthesis code is implemented into a Metropolis- 
Hastings Markov Chain Monte Carlo algorithm that eval- 
uates the posterior distribution function for the Bayesian 
inference. This posterior distribution results from taking 
into account the a-priori knowledge about the parame- 
ters and the information introduced by the observations. 
The names and abbreviations of the six parameters that 
describe the clumpy models, and the intervals consid- 
ered for the general fitting are detailed in Table [5] (see 
also Figured]). 

We introduced to the code the intervals of the param- 
eters reported in Table [6] as uniform distributions, thus, 
giving the same weight to all the values in each inter- 
val. The exception is the torus radial extent parameter 
Y, fixed as a Gaussian distribution centered at 15, with 
a width of 2.5. We fixed this parameter after confirm- 
ing that it is otherwise unconstrained in the fits and is 
not correlated with any of the other parameters^. The 
choice of F = 15 is justified below f ii4.3.ip . For those 
galaxies for which we have additional a-priori knowledge 
about any of the parameters (e.g., the inclination angle 
of the torus for galaxies with known masers), the priors 
have been set as narrow Gaussians centred at a given 
value. Apart from the six parameters that characterize 
a clumpy model, there is another additional parameter 
that accounts for the vertical scale displacement, which 
we allow to vary freely, needed to match the fiuxes of 
a chosen model to an observed SED. This vertical shift 
scales with the AGN bolometric luminosity (ij5)). 

In order to compare with the observations, 
BayesClumpy simulates the effect of the employed 
filters on the simulated SED by selecting the filter and 
introducing the observed flux and its corresponding error 
(assuming Gaussian errors). For a detailed description 
of the Bayesian inference applied to the clum py models 
see lAsensio Ramos fc Ramos Almeida! ()2009t ). 

4.3. Model Results 
4.3.1. Seyfert 2 Individual Fits 

The results of the fitting process are the marginal pos- 
terior distributions for the six free parameters that de- 
scribe the clumpy models and the vertical shift. These 
are the probability distribution of each parameter, rep- 
resented as histograms. Uniform priors have been em- 
ployed in our analysis of the observational data. There- 
fore, when the observed data introduce enough informa- 
tion into the fit, the resulting probability distributions 
will clearly differ from uniform distributions, either show- 
ing trends or being centered at certain values within the 
considered intervals. 

We fit the individual Sy2 SEDs with the 
BayesClumpy code, modelling only the torus emis- 
sion, assuming no direct AGN continuum is detected. 
We also assume that foreground extinction from the 
host galaxy does not significantly modify the infrared 
SEDs o f the Seyfert gala x ies in the sample. For ex- 
ample, iMartini fc Poggel ()1999D used V-H colors of 

^ The only exception is Centaurus A (which Y parameter results 
centered in ~10-15 without having constrained it. See Appendix 
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TABLE 6 

Clumpy Model Parameters and Considered Intervals 



Parameter 



Abbreviation 



Interval 



Radial extent of the torus 

Width of the angular distribution of clouds 

Number of clouds along the radial equatorial direction 

Power-law index of the radial density profile 

Inclination angle of the torus 

Optical depth per single cloud 



No 

1 
i 

■TV 



15 

15°, 75°' 
1, 15] 
0, 3] 
0°, 90°] 
10, 200] 



Sy2 galaxies in the CfA sample to derive levels of 
optical extinction up to Ay ~ 5 mag. Based on this, 
lAlonso-Herrero et alj (|2003t ) applied different values of 
A y up to 5 magnitudes to a subset of to rus models (those 
of lEfstathiou fc Rowan-Robinsonlll995f ) to illustrate the 
effects of dust in the host galaxy, and concluded that 
moderate amounts of foreground extinction {Ay < 5 
mag) will not have a significant effect on the nuclear 
fluxes at A > 2 /xm. Thus, we do not consider foreground 
extinction in the fits with the clumpy models for the 
Seyfert galaxies in our sample. The only exception 
is Centaurus A, since its core is h eavily obscured by 
a dust lane up to A y ~ 7-8 mag ( Hough et al.lll987l : 
iPackham etall Il996t iMarconi etahl |2000[) . Foi details 
on the specific fit of this galaxy see Appendix El 

The histograms corresponding to the probability dis- 
tributions of the free parameters are shown in Figure [5] 
for Circinus, which constitutes the best fit in our sam- 
ple. It also includes the histogram corresponding to the 
optical extinction produced by the torus along the LOS, 
computed as A^'^^ = 1.086 N( 



Ty e 



mag. 



We adopt a radial extent of F = 15 as a Gaussian 
prior with a width of 2.5 because this parameter is oth- 
erwise unconstrained and presents no correlation with 
other parameters. The chosen value of y=15 is consis- 
tent with all current observations of nearby active nuclei, 
which place the radial extent of t he torus within values 
iJaffe et al.l 12 004: Packham et al.l 
IVlcisenhci mer et al.l l20'07t 



-20-30 



small e r than 
20051: iTristram et al 



Radomski et al.ll200S : 



20071 



Raban et al.,, 20091) . and perhaps 



even as small as ~5-10. The lack of restriction of the 
torus radial ex t ent is in agreement with the findings of 
iNenkova et~aI1 (|2008bf ). who report that IR SED fitting 
sets only a poor constraint on the torus size (see Section 

ESI). 

From Figure [5] it is clear that the Circinus data (8 
photometric data points) provide sufficient information 
to constrain the model parameters. In addition to the 
above mentioned prior for the Y parameter, we introduce 
a Gaussian prior on the inclination angle of the torus i, a 
Gaussian centered in 85°, with a width of 2°. We make 
this restriction in the i paramete r because of the detec - 
tion of a water vapor megamaser (|Greenhill et al.ir2003D . 
which constrains the viewing angle to i 90°. The opti- 
cal depth per cloud results in a narrow Gaussian centred 
at the median value Ty = 30±2. The number of clouds 
results in a Gaussian-like distribution of median value 
A^o = 10±2, and the width of the angular distribution 
has a median value of a = 61°±8. Finally, we establish 
a lower limit of q > 2.4 at a 68% confidence level. The 
optical obscuration produced by the torus along the LOS 
would be of Af^^S = 320±fg mag. 



Although the solution to the Bayesian inference prob- 
lem are the probability distributions shown in Figure 
[5l we translate these results into corresponding spec- 
tra (Figurc[6]). The maximum-a-posteriori (MAP) values 
of the parameters (i.e., the modes) represent the "best 
fit" to the data, since the mode corresponds to the most 
probable value of each parameter, and we plot the corre- 
sponding spectrum as a solid line. The dashed line shows 
the spectral model obtained using the median value of 
the probability distribution of each parameter, which is 
characteristic of the observed SED. Finally, the shaded 
region indicates the range of models compatible with a 
68% confidence interval for each parameter around the 
median value. The observed Circinus SED is extremelly 
well-fitted with the clumpy models, including only the 
reprocessed torus emission without any stellar contribu- 
tion. The model fitting to the photometric points makes 
a spectroscopic prediction: the 10 fim silicate feature ap- 
pears in shallow absorp tion in the fitted models, as ex- 
pecte d for Sy2 galaxies (jllao et al.ll2007t iLevenson et al.l 
I2007D . In particular, the model prediction for Circinus is 
in qualitative agreement with spectroscop ic observations 
()Roche et al.l 120061 : ITristram et alll2007| i that show the 
silicate feature in absorption, although deeper than the 
predicted by the models (see Appendix E| • 

The excellent fit of the Circinus SED is the result of the 
combination of the good behaviour of the clumpy models 
and the accuracy of the IR nuclear fluxes. The near- 
diffraction- limited adaptive optics NACO/VLT ncar-IR 
data (providing resolution between 0.1"and 0.2") and the 
T-ReCS/Gemini South mid-IR data (resolution between 
0.3"and 0.5"), together with the proximity of the galaxy 
provide the least contaminated (by host galaxy emission) 
SED of our sample, and consequently the best example 
of an AGN-dominated SED. 

We apply the same technique to all the Sy2 individ- 
ual galaxy measurements (with the exception of NGC 
5728, for which we did not find near-IR nuclear fluxes in 
the literature, resulting then excluded from any further 
analysis), and plot the spectral results in Figure [71 The 
probability distributions of each parameter for all the in- 
dividual galaxies are summarized in Table [71 Comments 
on the individual fits of the Sy2 galaxies in our sample 
are reported in Appendix [Al Despite the more limited 
flux measurements in most individual cases, the resulting 
model fits are generally good. The fitted parameters are 
consistent with unified AGN schemes and the require- 
ment that the central engine be hidden from direct view 
in Type-2 AGNs. The only exceptional cases are NGC 
1808 and NGC 7582, as both exhibit intense nuclear star 
formation and have lower spatial resolution mid-IR mea- 
surements (obtained at the CTIO 4 m). Thus, we con- 
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Fig. 5. — Probability distributions of the free parameters that describe the clumpy models resulting from the fit of Circinus, together 
with the calculated value of the optical extinction along the LOS A^^^ . The vertical shift and the Y parameter have been marginalized. 
Solid lines represent the mode of each distribution, dashed lines correspond to the median, and dotted lines indicate the 68% confidence 
level for each parameter around the median. 



5F 



O 




10 

Wavelength [|j.m] 



100 



Fig. 6. — High spatial resolution Circinus SED. The solid line represents the "best fit" to the data, corresponding to the clumpy model 
described by the combination of parameters that maximizes the probability distributions of the parameters. The dashed line represents 
the model that corresponds to the median value of the probability distribution of each parameter. The shaded region indicates the range 
of models compatible with a 68% confidence interval for each parameter around the median. 



sider these two model results unreliable, especially that 
of NGC 7582. 

From the individual fits with the clumpy models, we 
infer the following conclusions. (1) For the majority of 
the Sy2 considered here, the clumpy models reproduce 
the observed SEDs, suggesting that the high spatial res- 
olution measurements are dominated by the reprocessed 
emission of the torus. (2) The 10 ^m silicate feature ap- 



pears in absorption in the fitted models for seven out of 
the eleven Sy2 (Figure [7|) and is practically flat in the 
case of NGC 7172. The unrehable NGC 1808 and NGC 
7582 results, and also the Centaurus A fit suggest silicate 
emission. This is probably related with a characteristic 
SED shape, characterized by low N/Q ratio plus L band 
excess (associated with stellar contamination). See Sec- 
tion 15.21 for details. (3) The average number of clouds 
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TABLE 7 

Parameters Derived from the Clumpy Model Fitting 



Galaxy Seyfcrt a Na q i n/ Ay^^ 

Type Median Mode Median Mode Median Mode Median Mode Median Mode Median Mode 



CGntSyiinis A 


2 


61it8 


62 


2±1 


1 


<0.2 


0.1 


>85 


89 




186 


300+iS'' 


260 


111 vc\ nut; 


2 


59 


10it2 


g 


>2.4 


3.0 




85 


30+2 


30 


55 

780+||5 
185+lf 


295 


IC 5063 


2 


>57 


75 


>11 


14 


<1.5 


0.4 


>65 


89 


70±|i 
30+™ 


66 


580 


Mrk 573 


2 


<37 


19 


6±2 


4 


<1.4 


0.1 


85 (fix) 


85 


27 


150 


NGC 1386 


2 


50±}i 


51 


11±3 


12 


1 c:+0.9 


2.1 


85 (fix) 


85 


47±g 


51 


<1460 


410 


NGC 1808 


2 


<35 


16 


8±t 


5 


1.2±0.9 


0.9 


37±i 


44 


198 


<140 


10 


NGC 3081 


2 


>52 


74 


10±3 


11 


>1.2 


2.9 


>42 


89 


34 


<450 


120 


NGC 3281 


2 


>68 


75 


9±l 
5±f 


5 


1 3+11 


0.6 


>49 


88 


<12 


10 


55+11 


50 


NGC 4388 


2 


>53 


74 


7 


<1.6 


0.2 


>38 


88 


<14 


10 


<80 


50 


NGC 7172 


2 


>54 


74 


5 


>1.7 


2.9 


>45 


89 


<12 


10 


50+20 


50 


NGC 7582 


2 


<29 


16 


<2 


1 


>2.5 


3.0 




58 


<27 


14 


<6 


1 


NGC 1365 


1.8 


35±^" 


31 


7±3 


3 


<1.7 


0.1 


<42 


7 


110+60 


88 


<170 


10 


NGC 2992 
NGC 5506 


1.9 
1.9 


45±i| 

25±y 


52 
15 


<2 


4 
1 


<1.0 

2 5+0-3 


0.1 
2.7 


>53 
85 (fix) 


85 
85 


36+if 
<68 


31 
22 


160+ig« 
<90 


130 
30 


NGC 3227 


1.5 


33±* 


23 


6±2 


3 


<2 


0.1 


<48 


14 


115+55 


77 


<210 


10 


NGC 4151 


1.5 


<32 


16 


<3 


1 




1.9 


4i±i 


52 


120+ti 


93 


<65 


10 



Note. — The Y parameter is fixed to be a Gaussian distribution centered at 15 with a width of 2.5. For the galaxies Circinus, Mrk 573, 
NGC 1386, and NGC 5506 the % parameter is also introduced as a Gaussian prior into the computations, centered at 85°with a width of 
2°, based on other observations. Probability distributions presenting a single tail have been characterized with the mode and upper/lower 
limits at 68% confidence. Syl.8 and Syl.9 have been fitted with the geometry corresponding to torus emission-only. The Syl.5 models 
include the intrinsic AGN continuum emission. The Sy2 NGC 5728 and the Syl NGC 1566 are not included here because of the lack of 
near-IR nuclear fluxes for them, preventing the SED fitting. 



along an equatorial ray is typically within the interval 
iVo = [5, 15]. The exceptions are Centaurus A (A^o = 
2±1) and NGC 7582 (A^o < 2). (4) High values of a are 
preferred from our fits {a = [50, 75]). We only find val- 
ues of cr < 40 for the galaxies Mrk 573, NGC 1808, and 
NGC 7582. (5) High values of the inclination angle of 
the torus are preferred, with i > 40°, in general. (6) The 
combination of z, ct, and A^o determines the likelihood of 
encountering a cloud along the line of sight to the central 
engine. In these cases, such encounters are always likely, 
which is consistent with having the AGN blocked from 
direct view, as expected in these Sy2. (7) Relatively low 
values of Ty (<100) are found. The exception here are 
Centaurus A and NGC 1808, for which we find Ty > 100. 
(8) The radial density profile q does not show any clear 
trend, and we find both low- and high-values within the 
considered interval. All the above mentioned intervals 
or limits of the parameters correspond to median values. 
The reason for choosing median values instead of modes 
is that the median gives a less biased information about 
the result, since it takes into account degeneracies, while 
the mode does not. 

4.3.2. Intermediate- Type Seyfert Individual Fits 

We fit the Syl.8 and Syl.9 nuclei the same way we fit 
the Sy2 SEDs, namely considering only the reprocessed 
torus emission. The optically broad lines and the rela- 
tively strong near-IR fluxes of the Syl.5, however, sug- 
gest that some direct AGN contribution is also present 
in these cases. Thus, wc include the intrinsic AGN emis- 
sion as a broken power law in these models. The AGN 
scales self-consistently with the torus flux, and additional 
extinction (separate from the clumpy torus) is a free pa- 
rameter. The resulting fits are plotted in Figure [51 The 
Syl NGC 1566 is not included in the following because 
of the lack of near-IR nucler fluxes in the literature for 
it. 

From the individual fits of the intermediate-type 



Seyfcrt galaxies with the clumpy models, we find sim- 
ilarities and differences with the Sy2 results. (1) Com- 
pared with the Sy2, lower values of a are preferred in 
the fits of both the Syl.8 and Syl.9 (ct = [25, 45]) and 
Syl.5 (ct < 35). (2) All the intermediate-type Seyferts 
result in a low number of clouds along the equatorial 
plane {Nq = [1,7]). (3) Low values of the inclination an- 
gle of the torus are found for Syl.5 (i < 50°) and also 
for the Syl.8 NGC 1365. (4) q does not show any clear 
trend for intermediate- type Seyferts. and both low- and 
high-values within the considered interval are found, as 
for Sy2. (5) Low values of the optical depth per cloud 
are found for the two Syl.9 galaxies (ry < 70), whilst 
for the Syl.8 and Syl.5, high values {tv >100) are pre- 
ferred. (6) In contrast with the wide range of values of 
the optical extinction produced by the torus for Sy2, all 
the intermediate-type Seyferts show Ay < 200 mag. (7) 
Including some obscured AGN contribution fits the near- 
IR excess of the observed SEDs of the Syl.5 NGC 4151 
and NGC 3227. The AGN is obscured by up to a few 
magnitudes at V in each case. (8) Except for the case of 
NGC 2992, for which both the data and the fit are very 
similar to those of Sy2 galaxies, the 10 /im silicate fea- 
ture appears either in weak emission or is absent in the 
intermediate-type Seyfert fits. As in Section 14.3.11 the 
values of the parameters above mentioned correspond to 
the median. Although the results from the modelling of 
Sy2 and intermediate-type Seyfert SEDs seem to point to 
different trends in their torus parameters, a larger set of 
objects is needed to clarify if these differences are really 
significant. 

The results of our SED modelling of Seyfert galax- 
ies are generally in g oo d agre ement with the find- 
ings of iNenkova et all ()2008bl ). who reproduce IR 
observations of both Type-1 and Typ e -2 Seyferts 
from the literature llSanders et al. Il989t I Elvis et al.l 



119941 : lAlonso-Herrero et al.l l2003l : IPrieto et al.l l2004l : 
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iMason et all l2006t iHao et all [2007t iNetzer et al.1 l2007f) 

using clumpy torus models. They find No ~ 5-15 clouds, 
Ty ^ 30-100, a ~ 30°-50°, and the angular distribution 
of the clouds having a soft edge (the only distribution 
we consider here), with q = 1 or 2. They reproduce the 
observational data with compact torus sizes {Y ^ 5-10), 
although they claim that the radial extent of the torus is 
poorly constrained from SED fitting. 

5. INTERPRETATION OF OBSERVATIONS AND MODEL 

FITS 

5.1. NIR Flux and Torus Inclination 

The intermediate-type Seyferts exhibit characteristi- 
cally higher H/N ratios, and the Syl.5 have flatter SEDs 
overall. Hot dust emission from the directly-illuminated 
faces of the clumps close to the central engine contributes 
to the near-IR emission in all these cases. In addition, 
the Syl.5 include a component of the direct AGN emis- 
sion, (i.e., the tail of the optical/ultraviolet power-law 
continuum), which strongly flattens their IR SEDs. 

In general, the relative near-IR flux depends sensitively 
on the torus inclination angle. The individual Sy2 show 
the same trend toward low near- to mid-IR ratios and 
therefore more edge-on views {i > 40° considering the 
median values of the i distributions, and i 90° con- 
sidering the modes). Only the unreliable measurements 
for NGC 7582 and NGC 1808 show high near- to mid-IR 
ratios and thus indicate more pole-on views. In the con- 
text of the clumpy models, the presence of a cloud along 
the line of sight, which may occur from any viewing an- 
gle, results in Sy2 classification. Cloud encounters are 
more likely at large inclination angle, and thus we gen- 
erally find larger values of the inclination angle for the 
Sy2. However, we do not suggest that all Sy2 are viewed 
exactly through the equatorial plane. This small sample 
includes a number of galaxies whose tori are known to be 
nearly edge-on (Circinus, Mrk 573, NGC 1386, and NGC 
5506). Hence there is a selection bias for highly inclined 
tori in our sample. 

The likelihood of a cloud encounter more completely 
depends on the combination of viewing angle, A^O: a-nd a. 
The preference for somewhat lower values of i, Nq, and 
a in the intermediate-type Seyferts suggests that they 
present fewer clouds along the line of sight, even allowing 
some direct detection of the AGN continuum in the case 
of the Syl.5. This configuration also increases the like- 
lihood of unimpeded views of some directly-illuminated 
cloud faces (i.e., those on the "back" side of the torus), 
which increases the near-IR flux. These selection effects 
may account for our flnding of l a rger va lues of a in the 
Sy2 sample than iNenkova et al.l ()2008bD report as typ- 
ical for all Seyfert types (~ 30°), based on population 
statistics of Seyfe rt 1 and 2 galaxies (jSchmitt et al.ll200ll : 
IHao et al.l[2005al) . 

5.2. Stellar Contamination in the Near-IR? 

One characteristic spectral shape the models cannot 
reproduce well exhibits both strong emission toward the 
long-wavelength end of the near-IR (around L band) and 
simultaneously a steep N/Q ratio. NGC 4388 is the most 
severe example of this problem; NGC 3281, NGC 5506, 
and NGC 7172 show similar effects. Because the for- 
mal uncertainties on the near-IR points are smaller than 
those at Q, the model fits tend to reproduce the ncar-IR 



measurements while undcrpredicting the emission at Q 
in these cases. (While we do not observe NGC 7172 at Q 
directly, the models would not accomodate the strong L 
band emission combined with a rise of the SED toward 
20 /im, assuming a typical N/Q ratio.) 

We can interpret these results either as a failure of the 
models or as evidence for near-IR contamination. If the 
former, the models require some modification to produce 
additional L band emission without otherwise altering 
the spectral shapes significantly. However, the fact that 
the most severe discrepancies arise in the galaxies that 
have the lowest resolution near-IR measurements sug- 
gests instead that stellar or other extended contributions 
can be significant even longwards of 2 fim. The L band 
observations of NGC 3281, NGC 4388, NGC 5506, and 
NGC 7172 were obtained with natural seeing on 3 and 4 
m telescopes, so even these "nuclear" fluxes cover scales 
up to 150 pc. 

The possibility of stellar or other contamination of the 
Q band fluxes cannot be discarded, especially those from 
the 4 m CTIO. Where strong stellar emission takes place 
in the nuclear region of the galaxies, the emission of 
hot dust heated by this intense star formati on is an im- 
portant contribution at ^20 /um (see e.g., iHelou et al.l 
120041: IMason et ani2007D . Alternatively, the NLR can 
contribut e more strong l y at l onger wavelengths, as the 
models of iGroves et al.l (|2006f l show. 

For the CTIO observations the resolution achieved in 
the Q band was ~ 1", which corresponds to physical 
scales up to - 100 pc (in NGC 1365, NGC 1808, and 
NGC 7582). On the other hand, for the Gemini data 
the resolutions are between 0.5" and 0.6" in the Q band. 
Depending on the distance of each source, this implies 
having resolution from ^ 10 pc (e.g., Centaurus A and 
Circinus) to ~ 150 pc (e.g., IC 5063, NGC 3281, and 
Mrk 573). Thus, the steep N/Q ratios that the models 
cannot correctly reproduce could be due to the combined 
effect of stellar contamination around both the L and Q 
bands, or NLR contamination that is preferentially more 
severe at Q. 

5.3. Torus Size 

While the IR SEDs do not constrain the size of the 
torus, they are consistent with the small torus of clumpy 
models, confined to scales less than 10 pc. Uniform den- 
sity models require the dusty torus to extend over large 
dimensions, t o provide cool dust that produces the IR 
emission (e.g.. lGranato fc Danesdfl99^ . Fundamentally, 
in smooth distributions, the dust temperature is a mono- 
tonic function of distance to the nucleus. In contrast, in a 
clumpy distribution, different dust temperatures can co- 
exist at the same distance , including cool dust at small 
radii (jNenkova et al.|[2002f ). so large tori are not neces- 
sary. Indeed, the flat probability distribution we found 
when allowing K to be a free parameter (and as small 
as 5) explicitly shows that small tori can produce the IR 
emission we observe. 

In the clumpy models, with a steep radial density dis- 
tribution {q = 2), the SED is never sensitive to the outer 
torus extent because the majority of the clouds are lo- 
cated very close to the nucleus for all values of Y. With 
flatter radial profiles, more clouds are located farther 
from the central engine. These model SEDs are then 
sensitive to y, but the variations are evident only at 
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wavelengths longer than 20 fim. Far-IR spectral obser- 
vatio ns would provide useful constraints (|Levenson et al.1 
|2007() . but high spatial resolution measurements will re- 
main critical because the far-IR emission from dust that 
stars heat increases. 

5.4. 10 fim Silicate Feature and Column Density 
Estimates 

The model fits to the photometric data yield spec- 
tral predictions. The 10 /xm silicate feature is of par- 
ticular interest. It is always weak (in emission or ab- 
sorption) in the clumpy models, unlike smooth-density 
distributions in which it can app e ar in deep absorption 
(iDullemond fc van Bemmell l2005l : iLevenson et all l2007l : 
iNenkova et al.l l2008bf ) . Fundamentally, the weak sili- 
cate feature arises in the clumpy model because both 
illuminated and dark cloud sides contribute to the ob- 
served spectrum. While most views of the Sy2 torus 
are through absorbing dark cloud faces, silicate emis- 
sion from some bright faces fills in the feature, making it 
shallower. The weak silicate absorption or emission ob- 
served in Seyfert spectra indicate that a clumpy medium 
is th e dominant source of obscuration (|Levenson et al.l 
|2007() . Indeed, at high spatial resolution, where AGN 
emission is practically isolated, the silicate absorption 
has never been observ ed to be deep in Seyfert galax- 
ies (jRoche et al.l |2006[) . The clumpy models repro- 
duce the general observed trends of silicate emission 
in quasars and absorption in obscure d Seyfert galax- 
ies (iRoche et al.l 119911: IWeedman etahl 120051: iHao et al.l 



2005bt ISie bcnmor gen et all l2005t ISturm et all l2005l : 

Levenson et al. 2Q03). They can also accomodate the 



exceptional o bservations of silica te absorption in unob- 
scured AGNs (iRoche et al.|[T99l and emission in buried 
AGNs (jMason et al.ll2009[ ). with different clump distri- 
butions. 

We quantify the model silicate feature with its 
strength, which we define as t'^qJ^jji = In(Fcont) — 
hi{Fcore), relative to the underlying continuum. We cal- 
culate Fcont using a spline fi t and Fcore a ll owing peak 
wavelength to vary, following [Si^ockv et al.l (|2008l) . The 
values of t^q^jj^ are reported in Table El The silicate 
feature appears in shallow absorption in the fitted spec- 
tra of seven out of the eleven Sy2. The remaining four 
sources are Centaurus A, NGC 7172, and the unreli- 
able cases of NGC 1808 and NGC 7582. All of them, 
except NGC 1808, have pubhshed mid-IR spectra that 
show the silicate feature in absorption, contrary to our 
results (see Appendix [K\ for references). Excluding NGC 
1808 and NGC 7582, the incorrect prediction of the 
silicate feature in very weak emission (practically fiat) 
for NGC 7172 is probably due to the stellar contami- 
nation of the seeing-limited near-IR data described in 
Section 15.21 Indeed, the lack of Q band data for this 
galaxy limits the effectiveness of the modelling. The 
case of Centaurus A is more complicated, since the 
nuclear fiuxes come from a < 10 pc region, but the 
clumpy models reproduce the silicate feature in emis- 
sion. This could be related to the controversial nature of 
the near- and mid-IR emission of this galaxy. According 
to several publications, synchrotron emission is signifi- 
cant or even dominant (|Bailev et al.lll986l: iTurner et al.l 
119921: IChiaberge erani2001l : iMeisenheimer et al.ll2007D T 



although iRadomski et al.l ()2008l ) argue instead that the 
synchrotron contribution to the mid-IR flux is small. For 
a detailed discussion of Centaurus A infrared emission 
nature and comparison with mid-IR spectroscopic data 
see Appendix 1X1 

Regarding the intermediate-type Seyferts, the silicate 
feature appears in weak emission or absent, except for the 
case of NGC 2992 (a Syl.9 with a SED very similar to 
those of Sy2. See Table EJ, for which it is in absorption. 
In general, all the predicted silicate features, either in 
emission or absorption, are weak, which is a consequence 
of the contribution of both illuminated and shadowed 
cloud faces that the clumpy models predict for all the 
geometries considered. 

For the OHMc dust extinction profile (|Ossenkopf et al.l 
[1991) . we can derive the apparent optical extinction us- 
ing A'^^ = T^Q^j^x 23.6 for the galaxies with the silicate 

feature in absorption (TableEl column 2). This value can 
be compared with the optical extinction along the LOS 
derived from our fits using the mode values of Ng, a, 
TV, and i parameters {Ay'~'^ , Table |7l). Modes are used 
here instead of medians in order to avoid lower limits 
that would make more complicated the comparison with 
AyPP. As expected, we obtain A^'^'^ » (see Table 
m) since on the basis of the clumpy models, the true opti- 
cal depth along LOS does not coincide with the apparent 
optical depth measured from the absorption silicate fea- 
ture (Ay^^). The directly illuminated faces of the clouds 
seen along the LOS refill the absortion silicate band (pro- 
duced by the obscuring material), thus producing lower 
apparent optical depth values. On the other hand, A^*^"^ 
takes into account the full model cloud distribution along 
the LOS. 

The columns of material implied in the X-ray absorp- 
tion (Table [H column 5) should be in principle compa- 
rable to or larger than those inferred from the fit of IR 
data with the models (Table |8l column 4). Thus, we esti- 
mated N^^^ from the optical depth along the LOS, using 
the r elationship Nk'^^ /A ^/^^ = 1.9x10^^ cm~'^ mag~^ 
from lBohlin et all ( 19781 ). Larger absorbing column den- 



sities derived from X-ray measurements {Nj^ ) com- 
pared with N^'-'^ suggest that the X-ray emitting regions 
are affected by a larger amount of extin ction than the 
mid-IR ones. As it was firstly noted by iGranato et al.l 
(|1997| ). dust-free absorption unrelated to the torus could 
be the main source of X-ray obscuration. Thus, the X- 
ray absorber would be gas located inside the dust sub- 
limation radius. This is the case for most of the sam- 
ple, with Nfj-^'^y /Njf^ ranging from 2 to 25 (exclud- 
ing the unrealistic NGC 7582 result). On the contrary, 
comparable values of Njf^ and (those of Cen- 

taurus A, NGC 7172, and NGC 5506) mean that the 
dust-free absorption in these galaxies is low. We derive 
unrealistic results for IC 5063 and NGC 2992, for which 
j^x-ray ^ ^LOS ^ rpj^.^ ^^^^j ^j^^ kuo wu in- 

trinsic X-ray variability of these galaxie s ( PoUctt a et al 
1996t ITurner et al.lll997l: lYaaoob et al.l [2007,: .GiUi et al 



20001 ). In the case of variable sources, we cannot directly 



compare the column density predictions from IR data 
modelling and X-ray measurements, since the data have 
not been taken simultaneously. 
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TABLE 8 
Silicate Feature Predictions 



Galaxy 


app 


fmas) 


y^LOS 

fmae) 


j^LOS 
(cm.~^) 


j^X-rays 

{'cm^2) 


Reference(s) 


Centaurus A 


-0.67 




260 


4.9x10^3 


1.5X1023 


a 


Circinus 


0.86 


20 


295 


5.6x10^3 


4.0x1024 


b 


IC5063 


0.94 


22 


580 


LlxlO^"* 


2.6xl023 


c 


Mrk 573 


0.44 


10 


150 


2.8x10^3 


>1.0X1024 


d 


NGC 1386 


0.41 


10 


410 


7.8x10^3 


1.5X1024 


e 


NGC 1808 


-0.21 




10 


1.9x10^2 


3.1Xl022 


f 


NGC 3081 


0.91 


21 


120 


2.3x10^3 


6.3xl023 


d 


NGC 3281 


0.21 


5 


50 


9.5x10^2 


2.0x1024 


g 


NGC 4388 


0.51 


12 


50 


9.5x10^2 


2.6xl023 


h 


NGC 7172 






50 


9. 5x10^2 


8.3xl022 


i 


NGC 7582 


-0.61 




1 


1.9x1021 


1.4X1023 


j 


NGC 1365 


0.00 




10 


1.9x10^^ 


4.8x10^3 


k 


NGC 2992 


0.74 


17 


130 


2.5x1023 


8.0x1021 


I 


NGC 5506 


-0.21 




30 


5.7xl022 


3.6xl022 


m 


NGC 3227 
NGC 4151 


-0.09 
-0.31 




10 

10 


1.9x10^^ 
1.9xl022 


5.3x10^^ 
6.9xl022 


n 
o 



References. — (a) [ Markowitz ct al. (2007); (b) Soldi ct al. (200 5i'l: (cllTurnc r ct al. (1997 ); (d) Th i s wor k; Ce ) ILevenson et al. (2006); 
CfllJimenez-Bailon etal (2005): (g) Vigriali &: Comastri ( 2002); Oij EIvis et~aTr(l2004) : (i ) Awaki et al.l (l200l): (i)[Tumer et al.l (2000): (k) 
IRisaliti et al.l I I2005I) : (D IYagoob et al.l ( feoOTi) : fm) lLrmer et al.l (feoOOl ): (n) ILamer et all (120031) : (o) IBeckmann et al.l (120051) ! 

Note. — Silicate strength t^q^^ measured from the fitted models; negative values indicate a silicate feature in emission. It was not 
possible to properly measure tJq'^jjj for NGC 7172 because the feature is in cxtrcmcUy weak emission, and self-absorbed. Correspond- 
ing apparent optical extinction, A'^^, assumes A'^^/t^q^^j^ = 23.6, and Ay'^^ comes from the modelling. N^'-'^ is calculated from 
Nh/A^'^^ = 1.9 X 102icm-2 mag-i (IBohlin et aU\W7^ . TV^"*""^ values are taken from the literature. 



5.5. AGN Luminosities 

The clumpy model fits yield the bolometric luminos- 
ity of the intrinsic AGN in each case. Combining this 
value with the torus luminosity, we derive the reprocess- 
ing efficiency of the torus. The vertical shift parameter 
scales with AGN luminosity, which we report in Table O 
(L^^^; column 1) for the median models. Integrating 
the corresponding model torus emission yields the torus 
luminosity, (L^°[; column 2). 

The Sy2 tori in our sample are efficient rcproccssers, 
absorbing and rc-cmitting the majority of the intrin- 
sic AGN luminosity in the IR (with mean efficiency 
0.6 ± 0.3). Among the Sy2, only Centaurus A, Mrk 573, 
and the unreUable results of NGC 1808 and NGC 7582 
indicate efficiencies below 50% (the mean efficiency is 
0.7 ±0.3, excluding NGC 1808 and NGC 7582). With 
the limited sample size, the efficiencies of the intermedi- 
ate types (0.4 ± 0.3) arc formally comparable to those of 
Sy2, but the measurements point out that possibly Type- 
2 tori reprocess radiation more efficiently than Type- 
1 tori, although the measured differences between both 
Seyfert types are not significant. A larger subset of Syl 
and intermediate-type Seyferts is needed to clarify if tori 
of different Seyfert types are genuinely different. The re- 
processing fraction depends primarily on the total num- 
ber of clouds available to absorb the incident radiation, 
which increases with the model parameters Nq and a. 
If all Seyfert nuclei are identical, only the viewing an- 
gle would change with the classification, not the proper- 
ties of the torus itself. While these results are limited, 
they suggest instead that classification is biased. Specif- 
ically, when the torus covers a smaller fraction of the sky 
viewed from the central engine (which results in a lower 
efficiency), the galaxy is less likely to be classified as a 
Sy2 (because the AGN is incompletely blocked from more 



lines of sight). Such a selection effect would also extend 
to Syl; IR observations of these galaxies would provide 
useful tests of this suggestion against strong forms of 
AGN unification. 

The outer size of the torus scales with the AGN bolo- 
metric luminosity: Ro = YRd, so assuming a dust sub- 
hmation temperature of 1500 K, Ro = 6(L^„'f ^/lO*^)^-^ 
pc (with fixed Y = 15). We find that all tori in 
our sample have outer extents less than 5 pc (Table 
[9]) , in agreement with recent mid-IR direct imaging of 
nearby Seyferts (jPackham et al.l l2005l : iRadomski et al. 
2008[) and also interferomet r ic observations (iJaffe et al. 
2004 'Tristram et aU l2007l : iMeisenheimer et al.l l2007r 
Raban ct al, 2009). Specifically, the estimated outer 
radii for Circinus {Ro = 0.6 pc) and Centaurus A 
(Rq = 0.4 pc) are clo se t o the upper hmits reported by 
iTristram et all ()2007f ) and IMeisenheimer et all (|2007| ). of 
1 and 0.3 pc, respectively. The uncertainty of these ra- 
dial sizes is large, mostly because the unknown value of Y 
enters the calculation linearly, but the consistency with 
other measurements further supports the assumption of 
Y = 15. 

The bolometric luminosity of the intrinsic AGN de- 
rived from the fits (Table |9l column 1) can be directly 
compared with the bolometric luminosities derived from 
the absorption-corrected 2-10 keV luminosities compiled 
from the literature (Table |9l column 5). To get ^j^^ 
from the intrinsic 2- 10 keV luminosit ies we applied a con- 
version factor of 20 (|Elvis et al.lll994D . In general, we find 
comparable values of ^jcbH ^'^d L(^f^^ within the consid- 
ered errors for L^^f^. Exceptions are NGC 3081, NGC 
4388, and NGC7172, for which L^fj)" ~ 10 x L^^^ , and 
the unreliable NGC 1808, in the opposite direction. 

6. CONCLUSIONS 
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TABLE 9 
BoLOMETRic Luminosity Predictions 






Galaxy 


T AGN 
(erg s 1) 


(erg s 1 ) 


T AGN 
'^Xbol 
(erg s-i) 


T AGN 1 T AGN 
'^Xbol 1 ^bol 



Ccntaiims A 


o.o 


X 


lU 


1 fi 
l.D 


X 


in42 


n 9Q 




1 .o 


X 


in43 


9 


Circinus 


1 n 


X 


iU 


Q 1 


X 




n c 1 
U.oi 


n fin 
U.OU 


1 o 


X 


in43 
iU 


1 o 
L.Z 


oUbo 


2.4 


X 


1 n44 
10 


2.1 


X 


1 n44 
lU 


0.87 




1. / 


X 


1 n44 
10 


U.7 


Mrk 573 


4.3 


X 




6.6 


X 


1043 


0.15 


3.93 


4.4 


X 


1044 


1.0 


NGC 1386 


3.4 


X 


1042 


2.0 


X 


1042 


0.59 


0.35 


1.3 


X 


1043 


4 


NGC 1808 


6.8 


X 


1042 


3.1 


X 


1042 


0.45 


0.49 


2.2 


X 


1041 


0.03 


NGC 3081 


1.0 


X 


10« 


8.9 


X 


1042 


0.89 


0.60 


1.0 


X 


1044 


10 


NGC 3281 


7.4 


X 


10« 


5.3 


X 


1043 


0.72 


1.63 


3.0 


X 


1044 


4 


NGC 4388 


2.1 


X 


10« 


1.7 


X 


1043 


0.81 


0.87 


1.5 


X 


1044 


7 


NGC 7172 


8.5 


X 


10*2 


6.7 


X 


1042 


0.79 


0.55 


1.1 


X 


1044 


13 


NGC 7582 


1.1 


X 


1044 


1.1 


X 


1043 


0.10 


1.99 


9.8 


X 


1043 


0.9 


NGC 1365 


1.7 


X 


104« 


1.1 


X 


1043 


0.65 


0.78 


3.0 


X 


1043 


1.8 


NGC 2992 


2.8 


X 


1043 


1.5 


X 


1043 


0.53 


1.00 


3.0 


X 


1043 


1.1 


NGC 5506 


7.3 


X 


1044 


5.1 


X 


1043 


0.07 


5.13 


2.2 


X 


1044 


0.3 


NGC 3227 


2.0 


X 


1043 


1.1 


X 


1043 


0.55 


0.85 


3.8 


X 


1043 


1.9 


NGC 4151 


1.4 


X 


1044 


2.2 


X 


1043 


0.16 


2.24 


1.7 


X 


1044 


1.2 



Note. — Bolometric luminosities corresponding to AGN luminosity and integrated flux of the model emission (i.e., reprocessed 
luminosity). Luminosities arc good to a factor of 2. Columns 4 and 5 correspond to the fraction of energy emitted by the AGN that 
is reprocessed by the torus (L^j^^/i^*^^) and outer radius of the torus calculated using L^*^^. Absorption-corrected 2—10 koV X-ray 
luminosities are taken from the literature (references in Table |5J. ^x'bol derived from 20 xL^'''^. 



Wc report subarcsecond resolution mid-IR fluxes for 
eighteen Seyfert galaxies. These nuclear fluxes, in com- 
bination with published near-IR measurements at com- 
parable resolution are used to construct spectral energy 
distributions corresponding to torus emission. We an- 
alyze and fit these SEDs with the clu mpy dusty torus 
models of iNenkova et al.l (|200a l2008al [bh. The main re- 
sults are summarized in the following. 

• The high spatial resolution mid-IR nuclear fluxes 
reported in this work provide a spectral shape of 
the individual SEDs that is different from that of 
large aperture data SEDs (on scales of a few arc- 
seconds). 

• The shape of the average Sy2 SED, constructed us- 
ing only pure Sy2 galaxies and either diffraction- or 
near-diffraction-limited data, rises steeply towards 
the mid-IR, with an IR slope (from ~1 /im to 18 
^m) am = 3.1 ± 0.9. 

• The individual Sy2 SEDs are typically steep 
through the IR, with aju ranging from 1.8 to 3.8. 
Considering separately the near-IR (from ^ 1 to ~ 
9 /im) and the mid-IR (~ 10 to ~ 18 fim) slopes, 
we find aNiR ^ cxmir, with unir = 3.6 ± 0.8 and 
cmir = 2.0 ± 0.2 measured from the average Sy2 
SED. 

• The IR SEDs of the intermediate-type Seyferts in 
the sample are flatter {ajji — 2.0±0.4 for Syl.8 and 
1.9 and aiR = 1.6±0.3 for Syl.5) and present larger 
N/Q band ratios than those of Seyfert 2 galaxies. 
The near-IR excess that flattens the SED of all 
intermediate-type Seyferts is due to the contribu- 
tion of hot dust from the illuminated faces of the 
clouds, while direct AGN emission is also impor- 
tant in the case of Syl.5. 

• The clumpy models of INenkova etlHI (|2008ai rbl) 

successfully reproduce the SEDs of both the Seyfert 



2 and the intermediate-type Seyferts. The mod- 
els accommodate the range of ajR observed in the 
former, while providing optically thick obscuration 
along the line of sight. 

• The IR SEDs do not constrain the radial extent of 
the torus, Y, and this parameter remains uncorrc- 
lated with the other parameters. Because the outer 
torus contains the coolest material, high angu- 
lar resolution measurements at wavelengths longer 
than 15 /im are needed to reveal significant varia- 
tions in the torus size. 

• In the individual Sy2 fits, we generally find the 
number of clouds within the interval Nq — [5, 15], 
the width of the angular distribution of clouds 
within a = [50°, 75°], high values of the inclina- 
tion angle of the torus (i=[40°, 90°]), and optical 
depths per cloud ry < 100. The radial density 
profile q does not show any clear trend within the 
considered interval ([0, 3]). 

• We find small number of clouds in the individual 
intermediate-type Seyfert fits {Nq = [1,7]). Low 
values of <t arc preferred from the flts of Syl.8 and 
Syl.9 (ct = [25°, 50°]) and even lower for Syl.5 
(ct < 35°). We require direct (though extinguished) 
emission of the AGN to reproduce the near-IR ex- 
cess observed in the SEDs of Syl.5. 

• Views of Sy2 are more inclined than those of the 
Syl.5 galaxies. More importantly, the larger values 
of No and ct in the Sy2 models suggest that these 
central engines are blocked from direct view along 
more lines of sight, in contrast to the intermediate- 
type Seyferts, which present more clear views of 
the AGN. Due to the limited size of the analyzed 
sample, these differences are not significant. 

• Although the fits are based on photometric data, 
the models predict the behavior of the 10 /im sil- 
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icate spectral feature. Most of the individual Sy2 
galaxies show shallow silicate absorption. In all but 
one of the intermediate- type Seyferts, the feature 
is either in shallow emission or absent. 

• The columns of material responsible for the X-ray 
absorption are larger than those inferred from the 
model fits for most of the galaxies in the sample, 
which is consistent with X-ray absorbing gas lo- 
cated within the dust sublimation radius, whereas 
the mid-IR flux arises from an area farther from 



the accretion disc. 

• The combination of large number of clouds and 
high values of a results in tori that more efficiently 
reprocess the incident nuclear radiation than those 
with lower values of these parameters. 

• In the models, the outer radial extent of the torus 
scales with the AGN luminosity, and we find the 
tori to be confined to scales less than 5 pc. 



APPENDIX 
INDIVIDUAL OBJECTS 

Below, we comment on the individual fits of the galaxies in our sample, including Type-2 Seyferts and intermediate 
types. Unless the opposite is indicated, all the reported values from the fits correspond to or are calculated from the 
median values. 

Centaurus A. — Centaurus A (NGC 5128) is a nearby Sy2 nucleus, with its core heavily obscure d by a d ust lane. This 
makes i t visible only at wavele ngths longwards ~ 1 /im ( Sch reier et "anil998l : lMarconi et al.ll2000[) . Indeed. [Hough et all 
(|1987f ). [Packham et al.l (|1996( ). and lMarconi et all ( )2000 ) determined a value for the extinct ion caused by the dust lane 
of Ay ~ 7-8 mag. However, the total obscuration on our LOS to the core was estimated bv lMeisenheimer et al.l ()2007l ) 
to be Ay ~ 14 mag, including the torus material. Due to the existence of the dust lane, that is probably affecting our 
nuclear fluxes (specially the near-IR ones) , we exceptionally have taken this foreground extinctio n into account f or the 
fit of Centaurus A with the clumpy models (Figure [7|) , by fixing Ay = 8 mag and applying the ICalzetti et al.l ()2000l ) 
law to the fitted models. 

The clumpy models reproduce the observed photometry very accurately, but not the broad absorption silicat e 
feature detected in mid-IR spectroscopy and interferometry (jSiebenmorgen et "all l2004at iMeisenheimer et al.ll2007l ). 
This mism atch in the fitted model could b e due t o the synchrotron radi ation contribution to the infrar ed fluxes 
claimed bvlBai lcv et al.' ('1986(). iTurner et al.' ('I992fl. IChiaberge et al.l ()2001f ). and IMeisenheimer et al.l (|2007f ). Indeed, 
IMeisenheimer ct al . (2007) find from mid-IR intcrferometric observations that the contribution of synchrotron emission 
to the nuclear e missio n of Centaurus A is of ^80% at ~8 fjm and ~60% at ~13 /xm. H owever, the results of 
iRadomski et al.l ()2008f ). including the lack of variabihty at ~10 fxm. (see iRadomski et alll2008l and references therein) 
and zero polarization at ~1 mm (although low millimeter polarization may be associat ed with low-luminosi t y sour ces; 
iPackham et"aI1ll996[ ) are inconsistent with a synchrotron source as that suggested by IMeisenheimer et al.l (|2007t) . A 
more detailed analysis of the SED of this galaxy, e.g., including an optically thin power-law in the fits to account for 
the synchrotron radiation can be done. However, this kind of detailed analysis of an individual galaxy is beyond the 
scope of this work. 

The probability distributions derived from the fit of the Centaurus A SED with the clumpy models clearly constrain 
the number of clouds A*o, the optical depth of each cloud ry, and the width of the angular distribution a. Nq resembles 
a Gaussian distribution centered at the median value Nq — 2±1. The ry histogram shows an asymmetric shape, with a 
tail towards low values and of median ry = 176 at a 68% confidence level. High-cr values are more probable than others 
(median value of 58° ±^5), and the same for the inclination angle of the torus {i > 85°). On the contrary, low- values 
of q are preferred {q < 0.2). The radial thickness of the tor us Y has been introduced to the code as a Gaussian prior 
centered in 15 with a width of 2.5, as explained in Section 14.21 Indeed, the preliminary fit of Centaurus A without 
restrictions in any of the model parameters was the only one that constrained the radial thickness of the torus, with a 
median value of its probability distribution oiY ^ 15. The calculated median value of the optical extinction produced 
by the torus is A\f"^ = 300±^[;" mag for this galaxy. 

Circinus. — The Ci rcinus galaxy is a nearby Sy2 with a prominent cone-shaped region revealed in the [O III] images 
(|Wilson et al.ll2000t l. Both the accuracy of the 8 high spatial res olution flux measurem ents and the proximity of the 
galaxy make the resulting fit excellent. Water vapor megamasers (jGreenhill et al.ir2003f ) constrain the viewing angle to 
i ~ 90°. Thus, in addition to the Gaussian prior for the Y parameter, we introduce a prior on z, a Gaussian centered 
in 85°, with a width of 2°^°. The resulting fit constrains most model parameters. The optical depth per cloud results 
in a narrow Gaussian centered at the median value Ty = 30±2. The number of clouds results in a Gaussian-like 
distribution (median value of Nq = 10±2), and the width of the angular distribution has a median value of ti = 61°±8. 
Finally, we establish a lower lim it of q < 2.3 at a 68 % confidence level. The optical obscuration produced by the torus 
would be A^°^ = 320±|g mag. iRoche et all (j2006( ) found the core of Circinus compact and obscured by a substantial 
column of silicate dust from T-ReCS/Gemini-South spectroscopy. They obtained a value for the silicate absorption 
depth of Tio^ni = 1-6. Our fitted models reproduce the silicate feature in absorption (see Table [8]), in qualitative 
agreement with the observations, but shallower (rfj^'^j-Q = 0.86) than that reported bv iRoche et al.l ()2006D . We find 

that the absorption column density inferred from X-ray measurements {N'^~'^°'^^ = 4 x lO^'' cm~^) is larger than that 

We choose this value of i to allow for a narrow Gaussian distribution within the considered interval [0°,90°]. 
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Fig. 7. — High spatial resolution SEDs of the Sy2 galaxies. Solid lines correspond to the model described by the combination of 
parameters that maximizes their probability distributions. Dashed lines represent the model computed with the median value of the 
probability distribution of each parameter. Shaded regions indicate the range of models compatible with the 68% confidence interval for 
each parameter around the median. Green dots correspond to seeing-limited near-IR data. 



20 



Ramos Almeida ct al. 




The IR Nuclear Emission of Seyferts on Parsec Scales 



21 



responsible for the mid-IR absorption {N^'^^ ^ 6 x 10^^ cm ^). This is consistent with the X-ray obscuring region 
being closer to the central engine than the mid-IR emitting dust. 

IC 5063. — This Sy2 nucleus is hosted by a merger remnant classified as an elliptical galaxy (|Colina et al.lll99"ll ). An 
obscuring dust lane partially covers the eastern side of the galaxy. A peculiarity of this sourc e is its strong radi o 
luminosity, which is two orders of magnitude larger than the typical values for Seyfe rt galaxies (IColina et al.lll99ll) . 
From an N-band adquisition imag e of IC 5063 obtained with T-ReCS /Gemini-South. lYoung et al.l (|2007f ) tentativelv 



found the nucleus of this galaxy slightly resolved. However, from our T-ReCS Si2 and Qa band imaging the nucleus 
appears unresolved, and consequently, we subtracted the corresponding PSF at the 100% level to derive the nuclear 
fluxes. From the fit with the clumpy models, we establish lower limits for three of the parameters [a > 57°, Nq > 11, 
and i > 65°) and an upper limit for the index of the radial density profile (q < 1.5). The only Gaussian- like distribution 
resulting from the fit corresponds to ry (median value of 70±2|). is calculated to be TSOitm mag. The galaxy 

presents one of the highest L^^^ of the sample (2.4x10*^ erg s~^, see Table|9]), and also a high reprocessing efficiency 
(87%), consistent with the large values of A'o and a we find. X-ray observations with the ASCA satellite allow to 
infer an intrin sic 2-10 keV luminosity of ~10'*^ erg s""'^, that is within the range expected for Seyfert 1 galaxies 
([Turner et al.l[T997l) . The silicate feature appears in absorption in the fitted models, in q ualitative agreement with 
mid-IR spectroscopic observations obtained with T-ReCS/Gemini-South (jYoung et al.ll2007l ). However, the absorption 
band appears deeper in our models (t^q^^ = 0.94) than in the mid-IR spectrum (rio^m = 0.33). 

Mrk 513. — Mrk 5 7 3 is optically classified as a classical Sy2 nucleus (jTsvetanov fc WalshI Il992l ). and 
iRamos Almeida et al.l ()2008f ) recently re-classify it as an obscured NLSyl galaxy, based on near-IR spectroscopy and 
X-ray archival data. The X-ray observations of Mrk 573 do not straightforwardly reveal the intrinsic luminosity of 
this AGN because it suffers from Compton thick obscuration (jGuainazzi et al.ll2005D . The strong soft X -ray emission 
is likely photoionized line emission. It is not the AGN continuum in any case , so previous simplified fits of power law 
models over a broad energy range (e.g.. lGuainazzi et al.ll2005l : IShu et al.ll2007[ ) fail to accurately measure the AGN. We 
determine the intrinsic AGN luminosity from the XMM-Newton observation of Mrk 573 with the pn detector, which 
offers a high signal/ noise ratio. We use the reduced data and associated background and calibration files from XAssist 
pipeline processing (|Ptak fc Griffiths! |2003[ ). We fit only the data at energies greater than 3 keV, to avoid confusion 
from the separate soft emission sources. The observed continuum is flat and the Fe Ka line equivalent width (EW) 
is large, which are characteristic of a purely reflected (not direct) AGN continuum. We fix the photon index F = 
and find EW = 2.0 ± .7 keV. The results ar e similar if a typical AGN continuum slope is adopted (F = 1.9). We 
follow the technique of iLevenson et all (j2006[ ). using the Fe line luminosity and EW and to determine the intrinsic 
AGN luminosity, finding Ljf^^ = 2.2 x 10^3 ergs ^ in the 2-10 keV bandpass. We note that the fitted hard continuum 
directly corresponds to an intrinsic luminosity that is a factor of 100 smaller. 

The IR SED shape of Mrk 573 is practically the same as those of the rest of Sy2 galaxies in the sa mple, so w e 
consider it equivalent to the other Sy2. This galaxy presents a strong symmetric double radio lobe (Nagar et al.lll999h . 
suggesting that it contains an edge-on torus; the radio lobes would be asymmetric if the jets were out of the plane 
of the sky. The SED fit reproduces all 7 data points and also constrains the majority of the parameters. As for the 
case of Circinus, we have introduced the inclination angle of the torus as a Gaussian prior, centered in 85° with a 
width of 2°. The probability distribution of TYq present a Gaussian shape, with median value iVo — 6±2 (there is a tail 
towards large values that enlarges the error) . The a and q distributions allow to establish upper limits of ct < 37° and 
q< 1.4. The optical depth of each cloud shows a Gaussian distribution with median value of ry = 30±g°. The optical 
extinction produced by this clumpy torus is Ay~^^ = 185±7§° mag. Mrk 573 has the highest L^^^ of the Sy2 subset 
(4.3 X lO^"' erg s~^; Table [5]), but due to the low values of A^o and <t derived from its fit, its reprocessing efficiency 
is low (^15%). The 10 /im sihcate feature appears in shallow absorption in the median fitted model {ti^J^jj^ = 0.44; 
see Table [U, producing an estimated value of the optical extinction A"^^ = 11 mag. The hydrogen column density 
obtained from the modelling along the LOS {N^'-'^ ~ 3 x 10^"^ cm~^) is much lower than the one that obscures the 
X-ray emission (A^^"™^'* > 1 x lO^^cm-^). 

NGC 1386. — NGC 1386 is one of the nearest Sy2 galaxies and therefore extensively studied. iMauder et al.l (jl992f ) 
resolved its NLR on a 0.3" scale (~ 15 pc) using speckle interferometry, detecting i ndividual NLR clou ds and claiming 
a clumpy structure. The detection of a water vapor megamaser in NGC 1386 (jBraatz et al.l fl997[ l constrains the 
inclination angle of the torus to be i ~ 90°, as in the case of Circinus and Mrk 573. The fit with the clumpy models 
reproduces the observed SED; we find A^o = H ± 3, cr = 50°±i[g, Ty = 95±5]', and q = 1.5±]';g at a 68% confidence 
level. The optical extinction produced by the clumpy torus is Ay'^^ < 1460 mag. The fitted models show a shallow 
silicate absorption feature (rfp^j^j = 0.41) from which we derived an optical extinction of A"^^ = 11 mag. As for most 

of the galaxies in the sample, we find N^'-'^ < N'^~^°'^^ ^ indicating that the X-ray obscuring region must be closer to 
the central engine than the torus material. 

NGC 1808. — The Sy2 galaxy NGC 1808 is undergoing an intense episode of star formation in its central 750 pc 
(|Tacconi-Garman et al.lll996l l2005l : iGaUiano fc AUoinl l2008f ) . The near-IR fluxes were obtained using ISAAC/VLT, 
and we measured the mid-IR flux from 10.7 and 18.2 /im OSCIR/CTIO images. Consequently, NGC 1808 has relatively 
low spatial resolution for our sample (the angular resolution of the ISAAC images is ^ 0.6" and that of the OSCIR 
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image at CTIO is ~ 0.9 — 1"). The near-IR SED shape is flat and peculiar, and star formation may contaminate these 
measurements (according to the argument of §5.2p . making the fit not very reliable. Low values of tr and i are preferred 
(< 35° and 37°±23, respectively), contrary to the general trend for the Sy2 galaxies in our sample. The median value 
of the number of clouds is A^o — Sif, the optical depth per cloud ry > 122, and q — 1.2±0.9. The calculated optical 
extinction produced by the torus is Ay"^^ < 140 mag. From the bolometric luminosites derived from the fit, we derive 
the reprocessing efficiency of the modelled torus, that results on ~45%. The 10 fj,m silicate feature appears in weak 
emission in the fitted spectrum (t°^^j^ — —0.21). Summarizing, due to the intense stellar formation that is taking 
place in the nuclear region of this galaxy, and to the lower spatial resolution near- and mid-IR measurements that we 
have for it (all from 3 or 4 m telescopes), the infrared SED of NGC 1808 is very likely contaminated with starlight. 
This is probably producing its peculiar shape and the different trends in the fitted parameters, in comparison with the 
rest of Sy2 in the sample. For these reasons, we consider the NGC 1808 SED fit unreliable. 

NGC 3081.—- NG C 3081 is a Sy2 galaxy with intense stellar formation located in a series of nested ringlike features 
()Buta et al.ll2004f ). As in the case of NGC 1386, the SED of NGC 3081 is formed with only three data points, and 
the model results are similar, with the probability distributions of the parameters presenting the same shapes. The 
A'o and ry histograms have median values of A'q ~ 10±3 and ry = 47±||, respectively. We determine lower limits of 
a > 52°, q > 1.2, and i > 42°. Ay^^ is estimated to be lower than 450 mag. According to the calculated luminosities 
for this galaxy, the fraction of reprocessed radiation (l^l"] / L-^^^) is among highest in the sample (~90%). The 10 
fim silicate feature is reproduced in absorption by the models (riO/Ltm°^^ = 0.91), from which we inferred an optical 
obscuration A"^^ ^ 23 mag. 

NGC 3281.— This Sy2 galax y suffers large extinction (Ay = 22 ±11 m ag to the IR-emitting and Nh = 2x lO^^'cm-^ 
to the X-ray emitting regions: [Simpsonlll998l : IVignali fc Comastrill2002l ). The clumpy models fit the data well, except 
for the 18.3 fiui point, which the models underestimate slightly. The probability distributions resulting from the fit 
show that large values of a and i are more probable than others (> 68° and > 49°, respectively) and establish an 
upper limit for the optical depth per cloud of ry < 12 at a 68% confidence level. The number of clumps Nq presents a 
Gaussian distribution with a tail towards large values, with median value A^o = 6±|, as the q histogram (median value 
of (? = 1.3±J g). The optical extinction produced by the clumpy torus corresponds to a Gaussian distribution of median 
value Ay'-'^ = 551b}! mag. The models predict the 10 fxm silicate feature in shallow absorption (riO/xm^P^ = 0.21), 
resulting in an apparent optical obscuration of A'^^ = 5 mag. Njf^^ < N'^~^°'^'^ , indicating a closer X-ray obscuring 
region to the central engine than the torus material. 

NGC 438 8. — This hig h ly inc lined disk galaxy is spectroscopically classified as a Sy2 nucleus (jPhillips et al.lfl98^ . 
However, [Shields et al.l ljl996[ ) report the detection of weak, broad Ha emission. The ground-based L-band point 
introduces a bump in the SED around 3-4 /itm, which likely represents starlight contamination (i j5.2p . The galaxy appear 
very extended in the L band (|Alonso-Herrero et al.l[2003h . making difhcult to isolate of the nuclear emission. Thus, 
the fit correctly reproduces the near-IR and the N-band points but severely underestimates the Q-band measurement. 
The probability distributions establish lower limits of cr > 53° and i > 38°, and upper limits oi q < 1.6 and ry < 14 
at a 68% confidence level. The median value of the number of clouds is A'o = 9±|. Along the line of sight. Ay < 80 
mag. The silicate feature appears in shallow absorption in the fitted models (riO/im°^'P = 0.51), allowing to estimate 
an optical obscuration of A"^^ = 13 mag. As for most of the galaxies in the sample, we found N^'-'^ < Njj~^'^^'^. 

NGC 7172. — NGC 7172 is a Type-2 Seyfcrt nucleus in a nearly edge-on spiral galaxy. A prominent dust lane in the 
east- western direction crosses the galaxy. Similar to the case of NGC 4388, we observe a near-IR bump of emission 
in the SED of NGC 7172, which stellar contamination of the seeing-limited fluxes likely produces ('i )5.2p . The Aq 
histogram shows a Gaussian shape, with median value Nq = High values of i, cr, and q arc more probable 

(i > 45°, a > 54°, and q > 1.7, respectively), while low values of ry are preferred (ry < 12). The torus produces 
optical extinction AyP^ = 50±20 mag. The reprocessing efficiency of this galaxy is high (~80%), according to the 
derived values of th e AGN and torus b olometric luminosities. T-ReCS mid-IR spectroscopy shows a substantial silicate 
absorption feature (jRoche et al.ll2007l ) that is not reproduced by the fitted models, that predict it in very extremely 
emission (and self-absorbed). This could be due to the above mentioned stellar contamination of the near-IR fluxes, 
that in addition to the lack of Q-band data point for this galaxy, make it difficult the correct reproduction of the 
silicate feature in absorption. 

NGC 7582. — NGC 7582 is a Sy2 galaxy that harbors a nuclear starburst surround ing the active nucleus 
(ICid Fernandes et al.l[200l ISosa-Brito et al.|[200lt IWold fc Gamanair2006l : iBianchi et al.ll2007t l and presents a nuclear 
outflow (iMorris et allll98 5n. The fit with the clumpy models shown in Figure [7] does not reproduce the 18.2 /xm point. 



and the resulting probability distributions of the parameters are very different compared with the rest of Sy2 galaxies 
discussed here (except for NGC 1808). Low values of a are more probable (< 29°), as are extremely small numbers of 
clouds (A^o < 2). High values of q (> 2.5) and low values of the ry (< 27) are preferred. The i probability distribution 
has a Gaussian shape which median value is i = 41°±2|. The optical obscuration produced by the torus results in 
an unrealistic value of Ay-'^ < 6 mag. The reprocessing efficiency fraction of this galaxy is among the lowest in the 
sample (^10%), a consequence of the low values of A^o and cr. The 10 /im silicate feature app ears in emission in 
the fitted models, in contradiction with TIMMI2 mid-IR observations ijSiebenmorgen et al.ll2004a[ ) that show a strong 
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absorption band characterized by rio^m = l-l- AH these results are probably due to stellar contamination of our 
mid-IR data (obtained at the 4 m CTIO). The intense circumnuclear star formation present in this galaxy makes it 
more difficult to isolate the torus emission from that of the host galaxy. Consequently, ss for the case of NGC 1808, 
we consider the model fit of NGC 7582 unreliable. 

NGC 1365. — This Syl.8 nucleus is hosted in a galaxy with intense nuclear star formation. Th e X-ray emission of 
NGC 1365 is also remarkable, showing the most dramatic spectral changes observed in an AGN (|Risaliti et al.l[2005l . 
[2001 . The rapid X -ray variability is attributable to variations in the LOS density, which is compatible with a low 
number of clouds along the LOS. The clumpy models fit the SED of NGC 1365 well, but the model parameters are 
not well-constrained. The number of clouds has a median value of 7±4 (which is consistent with the rapid X-ray 
variability), a = 35°±^§, and Ty = lll±|g. Low values of q and i (q < 1.7 and i < 42°) are more probable. The 
optical obscuration due to the torus is Ay'~^^ < 170 mag. The silicate feature is completely absent in the fitted models. 

NGC 2992. — This galaxy is optically classified as a Syl.9 nucleus, based on early publishe d spectra (IWard et al.l 



|1980[ ). which showed a weak broad Ha component, but no broad H/3. In a later observat ion (lAUen et al.lll999f) . the 
broad component of Ha dissapeared, leading to a new classification of the galaxy as a Sy2. iGilh et all (120000 reported 



that the spectrum regained the broad wings of Ha coincidin g with a period of i ntense X-ray activity, contrary to the 
low-stage that was taking place during Allen's observations. iTrippe et al.l (j2008[ ) report re-classification of NGC 2992 
as a Sy2 nucleus again, as of 2006. The near-IR ground-based data were taken in 1998 April, coinciding with the 
epoch of Gilli's observations, and thus with the classification as intermediate-type Seyfert. On the other hand, our 
mid-IR measurements from Michelle were obtained in 2006 May, coinciding with the Sy2 state. Thus, our IR SED 
includes fluxes taken during diferent periods of activity of NGC 2992. This can explain its similarity with a typical 
Sy2 SED. The clumpy models that include only torus emission (and no direct AGN component) fit the observed data 
well. The probability distributions are Gaussians for a, Nq, and ry, with median values a — 45° ±21, = 7±3, and 
Ty = 36±}f . High values of i (> 53°) and low values of q (< 1.0) are preferred. The calculated optical extinction 
due to the torus is Ay < 160 mag. NGC 2992 is the only intermediate-type Seyfert with the 10 /xm silicate feature in 
absorption in the fitted models (t^q^jj^ — 0.74). The derived apparent optical extinction from the measurement of the 

band is Ay'^ ~ 19 mag. The hydrogen column density along the LOS results in N^'^^ = 2.5 x 10^'^ cm~'^. This value 

results unrealistic when compared with that inferred from X-ray me asurements {N'^~'^ °'^^ = 8 x 10^^ cm~^). This 
could be due to the intrinsic X-ray variabihty observed in NGC 2992 (jTrippe et al.ll2008[ ). 

NGC 5506. — The nu cleus of this galaxy is classified as a Syl.9 b ased on the detect ion of broad wings of the Pa/3 
profile in the near- IR (jBlanco et al.lll990[ l. However, more recently, iNagar et al.l ()2002[ ) presented evidence that NGC 
5506 is an obscured Narrow Line Syl, as a result of the finding of the permitted O IA1.129 /xm line, together with a 
broad pedestal of Pa/3 and rapid X-ray variability. The clumpy models reproduce the near-IR and the N-band data 
points well, whilst slightly underestimating the Q-band flux. The discrepancy between the spectral observations and 
the model fit along with the poor mid-IR fit together suggest either that additional emission contaminates the near-IR 
measurements or t hat the models ca nnot produce this observed SED shape f ij5.2p . Based on the detection of a water 
vapor megamaser (iRaluv et al.lll998f) . we have introduced the inclination angle to the code as a Gaussian-prior centered 
in 85° with a width of 2°. In the resulting probability distributions, we find cr = 25°±y, q = 2.5±q-4, and the average 
number of clouds along an equatorial ray to be less than 2. The optical depth per cloud Ty < 68, and the calculated 
torus optical extinction is Ay < 90 mag. NGC 5506 has the highest Lf^^ of the sample (7.3 x 10^"* erg s~^; Table 
[9]). As a consequence of the low values of A^o and a derived from its fit, the reprocessing effi ciency is low (7%). The 
galaxy nucleus is very compact in the mid-IR, with an apparent optical depth of t^^^^ 1.4 (|R,oche et al.lll99lL [20071 : 
ISiebenmorgen et al.ll2004"b[ ) , although NGC 5506 also shows variations in its silicate absorption depth on parsec scales 
([Roche et al.l l2007( ) . Unfortunately, the fitted models reproduce the silicate feature in emission (t^q^jj^ = —0.21), 
contrary to the spectroscopic observations. This galaxy has comparable X-ray and mid-IR the absorbing columns, 
indicating that the dust-free absorption in this galaxy is lower than for the Sy2 reported here. 

NGC 3227. — This active nucleus is usually classified as a Syl. 5. XMM-Newton observations show short term hard 
X-ray variability likely related to variation in the intrinsic continuum emission ([Gondoin et al., 2003). The clumpy 
models reproduce the SED well, including the direct though extinguished AGN emission (using the lCalzetti "eraI1[2000l 
law), as in the case of NGC 4151. The resulting histograms show broad Gaussian distributions of the cr, Nq, and Ty 
parameters, with median values of 33°±f2, 6±4, and 116±55, respectively. Low values of i (< 48°) are more probable 
than others, and the q parameter results unconstrained. The optical depth calculated using the fitted model parameters 
is Ay'-'^ < 210 mag. On the other hand, the extinction that obscures the direct AGN emission is Ay = 2.7±2'5 '^^S 
(non related with the torus). This value is interm ediate among those rep orted in the literatu re of Ay = 1.2-1.7 
mag (|Cohenlll983[ : iGonzalez-Delgado fc Perezlll997[ ) and Ay = 4.5-4.9 mag ([MundeU et al.lll995f l. The 10 /xm sihcate 
feature is practically absent from the fitted models (t^q^jj^ = —0.1). 

NGC 4-151. — The nucleus of this well-known Syl. 5 galaxy sh ows flux variability in a wide wavelength range, 
with timescales ranging from a few hours in the hard X-rays ([Yaqoob et al.l Il993f l to several months in the IR 
([Oknvanskii et aL[|1999D . IShapovalova et ahl (|2008[ ) reported changes in the spectral type of the nucleus from a Syl. 5 
to a Syl.8, coinciding with maximum and minimum activity states, respectively. The SED of NGC 4151 is one of the 
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flattest in our sample, with the ncar-IR excess likely due to direct AGN emission, which is absent in Sy2 galaxies. 
We include this power-law contribution and allow for some extinction to account for the observed reddening of the 
active nucleus. The fit correctly reproduces the observed SED, especially in the near-IR. The probability distributions 
establish a median value of q = 1.7±Q g, i = 41°±|g, and ry = 120±4g. Low values for the number of clouds along 
equatorial rays (A'o < 3) and for the width of the angular distribution (ct < 32°) are more probable. The optical 
extinction derived from the fitted model parameters is Ay^^ < 65 mag. On the other hand, the optical obscuration of 
the direct AGN emission (not produced by the torus) is another free parameter of the fit, and we find Ay = 2.5±j | 
mag. The 10 /im silicate feature appears in wea k emission (t^q^^ ^ = —0.31) in the fitted models, in qualitative 

agreement with Spitzer observations of this galaxy (jWeedman et al.ll2005( l. 
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